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CHAPTER  I 

INTRODUCTION 

Geotextiles  have  increasingly  become  an  integral  part  of  many 
civil  engineering  applications  in  recent  years.  Slope  reinforcement, 
retaining  walls,  road  construction,  and  drainage  nets  are  a  few  among 
the  many  applications  of  these  materials  (3,12,16,19).  One  area  which 
has  been  receiving  more  and  more  attention  is  the  reinforcement  of  road 
embankments  overlying  soft  subgrade  soils  (3,12,16,18,19).  Much  of  the 
literature  relating  to  the  effect  of  geotextiles  upon  the  performance 
of  these  embankments  is  conceptual  in  nature  and  have  not  been  fully 
supported  by  experimental  and  theoretical  studies  (16).  These  analyses 
have  primarily  focused  upon  an  embankment  overlying  one  particular 
soft/weak  subgrade  soil.  The  analysis  of  an  embankment  overlying  more 
than  one  subgrade  soil  has  generally  not  been  required,  but  has  direct 
application  to  certain  areas  of  the  country. 

1 .1  Statement  of  the  Problem 

The  effect  of  frozen  ground  on  engineering  considerations  must 
allow  deviations  to  conventional  design  practices.  Normally  one  of  two 
broad  principles  can  be  followed  based  on  whether  or  not  the  frozen 
foundation  soils  are  thaw-stable  or  thaw-unstable  (14).  When 
foundation  soils  are  thaw-stable,  conventional  design  and  construction 
methods  may  be  used  (14).  Where  foundations  soils  are  thaw-unstable, 
conventional  design  and  construction  methods  must  be  modified  to 
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attain  sufficient  functional  and  economic  use  of  the  facility. 

Previous  studies  by  the  Alaska  Department  of  Highways  of  roadways 
constructed  over  permafrost  have  shown  that  progressively  deeper 
thawing  generally  occurs  beneath  annually  snow  covered  roadway  side 
slopes  (9).  The  insulating  effect  of  snow-covered  side  slopes  prevents 
the  complete  refreeze  of  the  underlying  soils  during  the  winter  which 
results  in  a  progressive  propagation  of  unfrozen  ground  or  talik  to 
occur  annually  (22).  The  degradation  of  this  talik  can  be  further 
accentuated  by  drifted  or  plowed  snow  accumulating  on  the  road-  way 
side  slopes  (14).  The  resulting  decrease  in  shear  strength  and  loss  of 
bearing  capacity  results  in  consolidation  and  slope  settlement  in  the 
thaw-unstable  permafrost  or  talik  (9,22).  The  thaw  settlement  which 
results  due  to  these  effects  usually  show  up  as  lateral  cracking  on  the 
wearing  surface  of  the  roadway  (22).  (See  Figure  1.1) 


Figure  1.1  Typical  roadway  distress  from  thaw  beneath  slopes. 

Methods  such  as  the  use  of  a  peat  underlay,  placement  of 
insulation  layers  within  the  embankment,  toe  berms,  and  air  duct 
systems  have  been  used  to  try  and  prevent  and/or  maintain  the 
foundation  soils  in  a  frozen  state  (9,10,22).  These  methods  have 


aided  in  the  short  term  maintenance  of  roadway  integrity  but  have  not 
solved  the  problems  of  long  term  thaw  related  settlements. 

The  use  of  geotextiles  as  a  means  of  facilitating  the 
construction  of  embankments  on  soft/weak  foundations  has  gained  in 
popularity  (16).  The  mechanisms  by  which  the  fabric  affects  the 
behavior  and  performance  of  conventional  aggregate-fabric-subgrade 
systems  include  the  following: 

(1)  Separation  of  the  aggregate  and  subgrade  soil 

(2)  Provision  of  a  filter  medium  to  facilitate  drainage 

(3)  Confinement  and  reinforcement  of  the  aggregate  layer 

(4)  Alteration  of  the  failure  mechanism  in  the  subgrade  soil  (19) 
The  application  of  these  principles  in  permafrost  regions  take  on  a 
slightly  altered  meaning.  Separation  of  the  aggregate  and  3ubgrade 

soil  would  primarily  be  affected  over  the  talik  zones.  Although  frozen 
soil  may  exhibit  deformation  and  separation  upon  loading,  the  area  of 
greatest  concern  would  be  over  the  thaw-unstable  talik  zone.  The 
fabric  overlying  the  thaw  zone  would  prevent  the  aggregate  and  thawed 
subgrade  soil  from  intermixing,  which  would  tend  to  reduce  the  effective 
depth  and  load  distributing  capability  of  the  aggregate  fill  (19). 

The  presence  of  a  geotextile  fabric  may  induce  the  formation  of 
a  natural  filter  in  che  subgrade  soil.  This  process  is  initiated  with 
the  migration  of  some  soil  particles  through  the  fabric  as  a  result  of 
fluid  flow.  Creation  of  the  filter  depends  upon  the  soil  grain  size  ' 
distribution,  the  porosity,  opening  size  and  opening  size  distribution 
of  the  geotextile,  and  the  hydraulic  gradient  of  the  flow  (19).  The 
hydraulic  gradient  may  be  influenced  by  the  in-si tu  stress  state  of  the 
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fabric,  where  tensile  and  compressive  stresses  may  alter  the  pore  size 
distribution  in  a  non-woven  fabric  (19).  The  filter  induced  by  the 
fabric  typically  maintains  the  separation  of  the  aggregate  and  subgrade 
soil  while  permitting  the  free  flow  of  water  (19).  When  the  frozen 
soil  under  the  toe  of  the  embankment  thaws,  the  resulting  water  is  not 
able  to  drain  off  properly.  The  surrounding  permafrost  does  not  lend 
itself  as  an  adequate  drainage  boundary  so  the  water  will  tend  to  flow 
through  the  fabric  as  the  outer  edge  of  the  embankment  consolidates  due 
to  the  loss  of  bearing  capacity  in  the  soil.  This  will  result  In 
differential  settlement  in  the  embankment  causing  longitudinal  cracks 
in  the  pavement.  This  has  been  one  of  the  primary  failure  mechanisms 
of  roadway  embankments  constructed  in  these  areas. 

The  primary  function  of  the  geotextile  is  provide  tension  rein¬ 
forcement,  transmitted  from  the  adjacent  soil  by  soil-geotextile  shear 
stress.  The  reinforcement  tensile  forces  resulting  from  shear  load 
transfer  transmit  soil  loads  from  one  location  (unstable  zone)  in  a 
soil  mass  to  another  location  (stable  zone)  (3).  Granular  materials 
have  zero  shear  strength  when  no  confinement  is  provided.  When 
confined,  granular  materials  furnish  excellent  load  bearing  capacity. 
The  load  bearing  capacity  of  the  aggregate  in  an  embankment  is  limited 
by  the  shear  stress  which  develop  at  the  aggregate-subgrade  interface 
(19).  A  layer  of  fabric  at  this  location  can  restrain  interfacial 
aggregate  movement  thereby  increasing  the  interfacial  shear  strength 
and  corresponding  load  bearing  capactiy  of  the  aggregate  layers  (19). 
The  load  bearing  capacity  of  the  aggregate  is  reduced  differentially 
when  talik  zones  develop  at  the  toe  of  roadway  embankments.  Additional 
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shear  stress  induced  by  the  geotextile  between  the  aggregate  and  talik 
zone  provides  confinement  for  the  aggregate,  thereby  increasing  the 
load  bearing  capacity. 


1 .2  Objective 

The  primary  objective  of  this  research  paper  is  to  estimate  the 
required  geotextile  properties  for  reinforcement  of  the  Farmers  Loop 
Road  in  Fairbanks,  Alaska  so  that  surface  deformations  can  be  kept 
within  tolerable  levels.  Numerical  analyses  are  used  to  evaluate  the 
effect  of  tensile  strength,  elastic  modulae,  and  location  of  the 
geotextile  on  the  surface  deformations  of  the  road. 

1 .3  Purpose 

The  purpose  of  this  research  is  to  provide  a  technically  sound 
and  economically  feasible  alternative  solution  to  the  long  term 
problems  associated  with  thaw-unstable  foundation  soils  underlying 
roadway  pavements  in  permafrost  regions  such  as  Fairbanks,  Alaska.  The 
functional  life  of  a  roadway  constructed  on  permafrost  is  decreased  by 
the  varying  properties  inherently  associated  with  frozen  soil. 
Substantial  changes  in  bearing  capacity,  shear  strength,  and  volume 
contribute  to  decreasing  the  functional  use  of  the  pavement.  Methods, 
previously  described,  have  been  attempted  to  increase  the  long  term 
functional  life  of  a  roadway  but  have  been  successful  only  in  the  short 
run.  The  main  purpose  of  this  project  is  to  evaluate  the  properties  of 
a  geotextile  which  may  be  used  to  reduce  the  magnitude  of  embankment 
deformations  given  the  progressive  degradation  of  the  permafrost 
beneath  the  embankment  slope  areas. 
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1  J4  Research  Approach 

The  success  of  any  method  of  analysis  largely  depends  upon  how 
well  the  actual  conditions  may  be  modelled.  For  embankments 
constructed  over  permafrost  four  basic  failure  modes  may  be  defined: 

(1)  Slope  failure  within  the  embankment. 

(2)  Rotational  failure  involving  the  embankment  and 
the  foundation  soil. 

(3)  Bearing  capacity  failure  of  the  thaw-unstable  talik. 

(i;)  Differential  settlement  of  embankment  as  a  result  of 

the  consolidation  of  talik  zone  under  the  side  slopes. 

Two  computer  programs  are  used  to  develop  the  required 
geotextile  properties  to  evaluate  stress  distributions  for  each  of 
these  failure  modes.  STABGM:  A  computer  program  for  slope  stability 
analysis  with  circular  surfaces  and  geogtextile  reinforcement  (7)  is 
used  to  evaluate  the  stress  distributions  for  the  first  two  modes  of 
failure.  Finite  Element  Analysis  for  Dams  (FEADAM)  and  Soil  Structure 
Interaction  Program  (SSTIPN)  are  used  to  analyze  stress  distributions 
and  evaluate  geotextile  properties  for  the  third  and  fourth  modes  of 
failure. 

Much  of  the  numerical  analyses  work  to  date  has  been  applied  to 
embankments  overlying  soft  soils  of  one  type.  This  literature  will  be 
used  to  develop  a  logical  approach  to  the  unique  problems  associated 
with  roadway  embankments  constructed  on  permafrost. 


CHAPTER  II 


SLOPE  STABILITY  ANALYSIS 

* 

The  computer  program  for  slope  stability  analysis  with  circular 
slip  surfaces  and  geotextile  reinforcement  (STABGM)  may  be  used  to 
analyze  internal  and  external  stability  of  roadway  embankments 
constructed  on  permafrost.  Internal  stability  refers  to  the  stability 
of  the  embankment  slope  alone.  Thus,  the  circular  slip  analyses  to  be 
employed  are  restricted  to  failure  through  the  toe  of  the  embankment. 
External  stability  refers  to  the  stability  of  the  embankment  and 
foundation  soil  combined.  For  embankments  constructed  or-  soft  soil, 
this  is  the  most  likely  mechanism  of  failure  (3). 

2.1  STABGM 

The  computer  program  STABGM  is  a  version  of  the  computer  program 
STABR  modified  in  December  1 98^4 ,  by  J.  M.  Duncan  and  B.  K.  Low.  In 
addition  to  the  slope  stability  analysis  capabilities  provided  by 
STABR,  STABGM  can  be  used  to  analyze  slopes  containing  up  to  twenty 
layers  of  reinforcement,  with  the  forces  in  the  reinforcement  layers 
varying  along  their  lengths. 

The  program  calculates  the  factors  of  safety  for  specified 
circles,  or  s marches  for  the  circular  slip  surface  having  the  minimum 


factor  of  safety,  using  the  Bishop's  Modified  Method  or  Ordinary  Method 
of  slices.  Briefly,  stability  analysis  using  the  ordinary  method  of 
slices  requires  trial  of  a  large  number  of  assumed  failure  surfaces. 

The  failure  zone  is  divided  into  a  series  of  vertical  slices.  It  is 
assumed  that  each  slice  acts  independently  of  its  neighbor.  There  is 
no  shear  developed  between  them,  and  the  normal  pessure  on  each  side  of 
a  slice  produced  by  the  adjoining  slices  are  equal  (2).  A  series  of 
calculations  are  then  evaluated  to  determine  the  resisting  moment, 
related  to  soil  strength,  and  the  overturning  moment,  related  to  the 
weight  of  the  soil  mass.  A  factor  of  safety  is  then  evaluated  for  each 
trial  failure  surface.  That  which  has  the  smallest  factor  of  safety  is 
the  most  critical  surface,  i.e.,  the  one  on  which  failure  is  most 
likely  tc  occur  (20).  The  Bishop's  Modified  Method  is  a  more  refined 
solution.  The  effect  of  forces  on  the  3ides  of  each  slice  are  taken 
into  account.  Through  a  series  of  trial  and  error  calculations,  a 
factor  of  safety  is  found.  The  ordinary  method  of  slices  yields 
results  which  are  too  conservative  (4).  The  Bishop's  Modified  Method 
is  more  representative  of  actual  conditions  and  is  used  for  the  slope 
stability  analyses  on  the  Farmers  Loop  Road.  The  program  may  be  used 
for  either  total  or  effective  stress  analyses,  or  a  combination  of 
both,  and  with  or  without  seismic  forces. 

2.1.1  Program  Operation 

The  program  consists  of  a  main  program  (STABGM)  and  three 
subroutines  (BISHOP,  ROOT,  and  SOLU).  Tie  geometry  of  a  slope  is 
describee  in  an  X-Y  coordinate  system.  The  X  coordinate  (horizontal) 
increases  from  the  top  to  the  toe  of  the  slope  and  the  Y  coordinate 
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(vertical)  increases  downward,  as  shown  in  Figure  2.1.  The  program  is 
capab?!.e  of  handling  irregular  slope  profiles,  tension  cracks,  soil 
layers  with  different  properties  and  nonuniform  thickness,  complicated 
pore  pressure  patterns,  and  irregular  variations  of  undrained  strength 
with  depth.  .... (7) 

Up  to  20  vertical  sections  may  be  used  to  define  the  geometry  of 
the  slope.  The  first  and  the  last  vertical  sections,  as  well  as  the 
last  soil  layer  boundary,  should  be  well  beyond  the  extent  of  any 
possible  slip  circle.  The  strength  of  the  soil  in  a  layer  may  be 
specified  either  by  the  cohesion  intercept  (c)  and  the  friction  angle 
( <t>)  which  are  constant  throughout  the  layer,  or  by  a  curve  of  undrained 
it  l/I*  V  crsus  depth.  The  total  number  of  points  describing  the 

variation  of  undrained  strength  with  depth  may  not  exceed  20 . (7) 

If  a  search  for  the  critical  circle  is  desired,  the  program 
searches  for  che  circle  having  the  minimum  factor  of  safety  by  either 
the  Ordinary  Method  of  Slices  or  Bishop's  Modified  Method.  Either  a 
horizontal  line  to  which  all  circles  are  tangent  or  a  specific  point 
through  which  all  circles  pass  must  be  defined.  In  addition, 
coordinates  of  the  center  of  the  first  circle  to  be  analyzed,  and  the 
final  grid  spacing  desired  in  the  search  must  be  specified.  The  search 
starts  with  calculation  of  factors  ol  safety  for  the  specified  circle 
center  and  for  centers  spaced  symmetrically  around  the  specified 
center,  as  shown  in  Figure  2.2.  The  centers  are  generated  in  the  order 
shown,  by  rotating  around  the  specified  center  with  a  spacing  whose 
length  is  equal  to  twice  the  final  grid  spacing.  ....(7) 

If  a  factor  of  safety  less  than  that  at  the  center  of  rotation 


•  FIG.  2.2  FIRST  CLOCKWISE  ROTATION  AROUND  THE  GIVEN  CENTER  A. 

THE  ROTATION  STARTS  AT  POINT  1 ,  WITH  RADIUS  OF 
ROTATION  TWICE  THE  FINAL  GRID. 


FIG.  2.3  THE  45-DEGREE  CLOCKWISE  ROTATION  AROUND  THE  CENTER  B. 
THE  ROTATION  STARTS  AT  POINT  1 ,  WITH  RADIUS  OF 
ROTATION  1.414  TIMES  THE  FINAL  GRID. 


is  found  at  any  point,  this  point  becomes  the  new  center  of  rotation. 

If  the  difference  between  the  factors  of  safety  at  the  old  and  the  new 
centers  of  rotation  is  equal  to  or  less  than  5%  of  the  factor  of  safety 
at  the  new  center  of  rotation,  the  spacing  is  reduced  to  the  final  grid 
spacing.  If  a  full  rotation  is  completed  and  no  factor  of  safety 
smaller  than  that  at  the  center  of  rotation  has  been  found,  a  new 
clockwise  rotation  round  the  same  center  is  initiated  with  the  spacing 
being  reduced  to  the  final  grid  spacing.  ....(7) 

If  a  full  rotation  with  the  final  grid  spacing  is  completed  and 
no  factor  of  safety  smaller  than  that  at  the  center  of  rotation  has 
been  found,  then  another  clockwise  rotation  around  the  same  center  is 
initiated.  This  rotation,  however,  starts  at  an  angle  of  H5°  with  the 
horizontal.  The  increments  of  rotation  are  still  90°,  thus  filling  in 
the  corners  of  the  final  grid,  as  shown  in  Figure  2.3.  If  no  smaller 
factor  of  safety  is  found  in  this  rotation,  the  factor  of  safety  at  the 
center  of  rotation  is  assumed  to  be  the  minimum  factor  of  safety,  and 
the  search  is  terminated.  ....(7) 

If  the  minimum  factor  of  safety  has  not  been  found  after  51 
circles  have  been  analyzed,  the  search  is  discontinued  without  printing 
out  a  minimum  factor  of  safety.  The  program  then  proceeds  to  search 
for  the  minimum  factor  of  safety  for  the  next  depth  limiting  tangent. 
When  this  occurs,  the  user  should  rerun  the  problem,  specifying  a  new 
starting  center  for  the  search.  It  is  frequently  desirable  to  specify 
a  rather  large  grid  spacing  for  the  initial  search,  so  that  the  center 
of  the  critical  circle  can  be  located,  at  least  approximately,  by 
analyzing  relatively  few  circles.  Then  a  finer  grid  spacing  may  be 
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specified  for  a  final  search  to  locate  the  center  of  the  critical 
circle  with  a  higher  degree  of  accuracy.  ....(7) 

With  STABGM,  the  improvement  in  the  factor  of  safety  is 
considered  to  be  due  to  the  increased  resisting  moment  provided  by  the 
reinforcement: 

4  mr  •  j,  tv**!-  y)  (i> 


in  which  AMr  »  resisting  moment  due  to  reinforcement 

=  force  in  i  th  layer  of  reinforcement  where 
it  is  intersected  by  the  slip  circle  being 
analyzed 

Yf^  **  Y  -  coordinate  of  reinforcement  (down  is 
positive) 

Yc  *  Y  -  coordinate  of  circle  center 
N  =  number  of  layers  of  reinforcement 
The  factor  of  safety  of  the  reinforced  slope  is  calculated  as 
follows: 


(2) 


in  which  Fw  =  factor  of  safety  for  a  given  circle  with 
reinforcement 

Fwo  =  factor  of  safety  for  the  same  circle  without 
reinforcement,  as  calculated  by  the  Ordinary 
Method  of  slices  or  Bishop's  Modified  Method 


AMr  =  resisting  moment  due  to  reinforcement,  calculated 
using  Equation  i 

M0  =  overturning  moment  for  the  circle 
STABGM  will  search  automatically  for  the  circle  with  the  minimum 
value  of  Fw  to  locate  the  circle  with  the  minimum  factor  of  safety. 

Up  to  twenty  horizontal  layers  of  reinforcement  can  be 
soecified,  with  as  many  as  nine  pairs  of  values  of  X-coordinate  and 
force  for  each  reinforcement  layer  to  describe  the  variation  of  force 
along  its  length. 

2 . 2  Reinforcement  of  Embankments  on  Weak  Foundations 
Geotextile  reinforcement  may  be  use  at  the  interface  between 
embankments  and  weak  foundations  to  prevent  failure  of  the  foundation 
soil  (3).  Two  performance  criteria  usually  considered  for  normal 
embankments  constructed  on  weak  soil  are  adequate  stability  and 
acceptable  total  and  differential  settlement  (3).  The  failure 
mechanisms  associated  with  embankments  constructed  on  permafrost  are 
similar  in  principle  but  not  in  sequence.  Initially,  embankments 
constructed  on  permafrost  are  very  stable.  Permafrost  is  relatively 
strong  as  a  foundation  soil,  but  once  thawed,  the  soil  becomes  very 
weak  relative  to  its  former  strength.  Additionally,  thawing  usually 
occurs  under  the  side  slopes  of  the  embankment  and  not  under  the 
centerline  which  introduces  the  problem  of  nonuniform  foundation  soil 
properties  and  differential  settlement  after  construction. 

Internal  stability  can  be  impaired  if  the  embankment  itself 
lacks  internal  stability  or  if  embankment  failure  results  from 


excessive  deformation  or  failure  of  the  foundation  soil  under  the 
embankment  load  (3).  Height  of  the  embankment  and  the  steepness  of  the 
side  slopes  govern  the  internal  stability  of  an  embankment.  Both 
parameters  also  influence  the  magnitude  and  distribution  of  stresses  on 
the  foundation  soil  and,  therefore,  influence  the  stability  of  che 
foundation  soil  (3).  This  interrelation  is  directly  applicable  to 
embankments  constructed  over  permafrost  which  has  thawed.  Upon 
thawing,  the  talik  zone  under  the  side  slopes  of  the  embankment  loses 
bearing  capacity  resulting  in  settlement  of  the  overlying  side  slopes. 
Consequently,  the  fill  material  may  weaken  and  loosen  thereby  impairing 
embankment  internal  stability  (3).  Embankment  internal  stability  can 
be  improved  by  layers  of  reinforcement  placed  within  the  embankment 
(3,11,12,16,17,18,19). 

The  placement  of  a  layer  of  reinforcement  at  the  embankment/ 
foundation  soil  interface  affects  embankment  internal  stability  two 
ways.  Upon  thawing  of  the  frozen  soil  under  the  side  slopes,  the 
embankment  soil  may  slide  along  the  geotextile  .nterface.  This 
mechanism  is  more  likely  to  occur  if  the  inter-  ;ce  friction  angle  and 
adhesion  are  relatively  low  between  the  geoteyt  _.e  and  the  embankment 
soil  (3).  If  lateral  sliding  does  not  occur,  the  reinforcement  layer 
will  aid  in  reducing  lateral  deformations  of  the  embankment,  thereby 
minimizing  the  risk  of  creating  voids  in  the  fill  from  cracking  (3). 

In  conventional  design  methodology,  prevention  of  embankment  failure 
due  to  lateral  sliding  is  the  first  step  (3,12). 

The  next  step  is  to  address  foundation  soil  stability.  External 
failure,  or  failure  of  the  foundation  soil,  is  said  to  occur  as  a 


result  of  one  of  two  mechanisms:  slip  surface  failure  or  bearing 
capacity  failure  (3).  Slip  surface  failure  occurs  when  a  portion  of 
the  embankment/ foundation  slides  relative  to  the  other  stable  part.  If 
reinforced,  the  slip  circle  would  likely  go  through  the  embankment,  the 
reinforcement,  and  the  foundation  soil.  A  bearing  capacity  failure 
would  result  when  the  embankment  punches  through  the  foundation  soil. 
This  would  occur  in  conventional  cases  when  the  reinforcement  in  the 
embankment  is  strong  enough  to  hold  the  embankment  together  to  prevent 
a  slip  surface  from  developing  (3).  In  permafrost  regions,  the  dynamic 
effect  of  frozen  soil  thawing  under  an  existing  embankment  does  not 
lend  itself  to  such  a  structured  approach.  Failure  of  an  embankment  in 
these  regions  is  likely  to  encompass  a  combination  of  all  three  failure 
mechanisms:  laterax  sliding,  slip  surface  failure,  and  bearing 
capacity  failure. 


2.2  Initial  Analysis 

Verification  of  embankment  stability  as  3tated  earlier  must 
consider  three  failure  mechanisms.  For  embankments  constructed  on 


permafrost,  these  failure  mechanisms  are  interrelated.  Once  the  frozen 
soil  thaws  under  the  side  slopes  of  the  embankment,  bearing  capacity  is 
lost  within  the  thaw-unstable  talik.  As  the  talik  consolidates, 
failure  planes  may  develop  within  the  embankment.  If  a  geotextile  is 
present  at  the  embankment/subgrade  interface,  the  risk  of  lateral 
sliding  on  the  outer  part  of  the  embankment  increases  and  may  add  to 
the  development  of  failure  planes. 

The  evaluation  of  internal  and  external  stability  of  the 
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embankment  was  initially  evaluated  using  STABGM.  The  analyses  of  toe 
ana  deep  stability  were  evaluated  based  upon  none,  one,  two,  three,  and 
four  layers  of  reinforcement  placed  within  the  embankment.  Analysis  of 
the  embankment  with  no  layers  of  geotextile  reinforcement  was  also 
conducted  to  determine  at  what  depth  the  most  critical  slip  circle 
surface  would  be  located  (i.e.,  the  depth  of  slip  circle  at  which  the 
lowest  factor  of  safety  was  evaluated). 

In  addition  to  critical  depths  of  slip  circles,  three  distinct 
dimensions  of  30ft  thaw  zone  subgrade  soil  were  also  analyzed.  The 
geometry  of  the  soft  thawed  zone  which  created  the  most  critical 
conditions  (i.e.,  lowest  factor  of  safety)  in  conjunction  with  the 
critical  slip  circle  depth  was  then  used  throughout  the  remainder  of 
the  analysis.  Once  this  critical  depth  and  critical  dimension  of  thaw 
zone  was  established,  further  analysis  of  stability  with  geotextile 
reinforcement  was  conducted. 

Ideally,  reinforcement  layers  should  be  placed  in  the  direction 
of  maximum  tensile  strain  (3).  For  steep  slopes,  principal  compressive 
strains  are  nearly  vertical  and  principal  tensile  strains  are  nearly 
horizontal  (3).  Therefore,  reinforcement  layers  placed  horizontally 
provide  a  high  degree  of  efficiency.  Also,  reinforcement  is  generally 
placed  in  horizontal  layers  to  accommodate  actual  construction 
sequences. 

2.3.1  Assumptions 

Figure  2.4  depicts  the  conditions  of  the  initial  analysis.  The 
following  assumptions  represent  the  initial  analysis: 

1.  Height  of  embankment  =  5  feet. 
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2.  Slip  circles  analyzed  were  located  tangentially  at  depths 
of  3,  5,  and  10  feet  below  the  bottom  of  the  embankment. 

3.  Thaw  zones  (soft  soil)  were  located  10,  15,  and  20  feet 
from  centerline,  symmetric  about  the  toe  of  the  embankment. 

4.  No  seismic  activity  was  taken  into  account. 

5.  No  pore  pressures  were  taken  into  account  (Total  stress 
analysis) . 

6.  Reinforced  embankment  stability  analyses  were  conducted  with 
geotextile  layers  located  at  5,  4,  2,  and  0.5  feet  beneath 
the  top  of  the  embankment. 

7.  Material  properties  were  modelled  as  crushed  stone  in  the 
embankment,  soft  clay  in  the  thaw  zone,  and  stiff  clay  in 
the  permafrost  zone. 

The  geometry  of  the  talik  zones  used  in  the  analyses  is  con¬ 
sistent  with  measurements  obtained  by  the  Alaska  Department  of  Trans¬ 
portation  and  Public  Facilities  (9,15,22).  The  material  properties 
were  modelled  according  to  relative  strength  to  each  other.  They  do 
not  represent  actual  material  strengths.  The  embankment,  crushed  rock, 
was  modelled  as  a  complete  unit  and  not  separate  layers  of  fill.  The 
thaw  zone  was  modelled  as  material  with  very  little  strength  relative 
to  the  embankment  and  permafrost.  Material  properties  selected  are 
shown  in  Figure  2.4. 

Embankment  stability  was  initially  analyzed  with  no  layers  of 
geotextile  reinforcement.  This  was  to  evaluate  at  which  depth  and 
which  dimensions  of  slip  circles  and  thaw  zones,  respectively ,  provided 
critical  conditions. 
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Analysis  was  then  conducted  with  layers  of  geotextile 
reinforcement  within  the  embankment.  Embankment  stability  was  then 
analyzed  with  one  layer  of  reinforcement  placed  at  the  embankment/sub¬ 
grade  interphase.  Additional  analyses  with  two,  three,  and  four  layers 
of  geotextile  placed  at  2,  and  0.5  feet  beneath  the  embankment, 
respectively,  were  conducted  to  determine  the  impact  on  embankment 
stability.  The  geotextile  reinforcement  strength  for  the  analyses  with 
one,  two,  and  three  layers  was  5,000  lbs  per  foot  of  width  for  each 
layer.  The  strength  of  the  upper  geotextile  layer  in  the  analysis  with 
four  layers  of  reinforcement  was  estimated  to  be  1667  lbs  per  foot  of 
width.  Mobilization  of  the  full  force  in  the  geotextile  layer  close  to 
the  surface  is  not  practical  due  to  the  minimal  overburden  pressure  of 
the  soil.  Also,  mobilization  of  the  full  force  in  all  four  layers  was 
modelled  to  occur  three  feet  from  the  end  of  each  layer  ( i . e . ,  zero 
force  applied  at  endpoints  and  5,000  lbs  force  applied  at  three  feet  in 
from  each  end  for  layer  one). 

2.3.2  Results  of  Initial  Analysis 

The  analysis  with  no  layers  of  reinforcement  provided  two 
important  pieces  of  information  which  impacted  the  remaining  analyses. 
First,  the  depth  at  which  the  critical  slip  circle  occurred  was  10  feet 
below  the  embankment/subgrade  interphase  (i.e.,  lowest  factor  of 
safety)  (See  Appendix  A).  Second,  the  dimensions  of  soft  soil  (thaw 
zone)  at  which  critical  slip  circle  conditions  occurred  were  at  both  10 
and  15  feet  from  the  centerline,  symmetric  about  the  toe  (i.e.,  factors 
of  safety  were  equal  for  both)  (See  Appendix  A). 
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The  analysis  with  three  and  four  layers  of  reinforcement 
revealed  that  no  appreciable  gain  in  safety  occurred  as  a  result  of  the 
fourth  layer. 

2.3.3  Revisions 

The  following  revisions  were  made  to  the  initial  conditions 
based  upon  the  results  obtained  (See  Figure  2.5): 

1.  Embankment  depth  was  increased  to  seven  feet. 

2.  Analysis  was  conducted  with  the  thaw  zone  located  10  feet 
from  centerline,  symmetric  about  the  toe  of  the  embankment. 

3.  Analyses  were  conducted  with  none,  one,  two,  and  three 
layers  of  geotextile. 

4.  Stability  analyses  were  conducted  with  slip  circles  located 
at  a  depth  of  10  feet  below  the  embankment  (external 
stability)  and  through  the  toe  of  the  embankment  (internal 
stability) . 

The  embankment  depth  was  increased  to  seven  feet  to  allow  for  a 
more  conservative  analysis.  The  greater  the  embankment  height,  the 
lower  the  factor  of  safety.  Analyses  of  external  stability  were 
conducted  to  establish  the  required  geotextile  reinforcement  to  achieve 
a  minimum  factor  of  safety  of  1 .0.  Analyses  of  internal  stability  were 
only  conducted  to  compare  factors  of  safety  with  external  stability 
analyses. 


Results  of  STABGM  Analyses 


The  results  of  the  revised  analyses  are  summarized  in  Table  2.1. 


Geotextile  reinforcement  strengths  of  5,000,  7,000,  9,000,  and  8,000 
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Table  2.1.  Summary  of  Results 


:• 


TANGENT  SLIP  CIRCLE: 

• 

Number  of 

Reinforcing  Layers 

Slip  Circle 

Depth  Below 
Embankment 

Factor  of 
Safety 

Force  of 
Layer 

F 

None 

3  ft 

5  ft 

10  ft 

0.719 

0.636 

0.593  (critical)- 

• 

One 

(7'  below  surface) 

10  ft 

1 .043 

3000  lbs 

Two 

(4*  and  7'  below  surface) 

10  ft 

10  ft 

1.035 

1 .035 

5000  lbs 
(each) 

:• 

Three  10  ft 

(2',  4'  and  7'  below  surface) 

1 .062 

1666  lbs 
(Layer  3) 
5000  lbs 
(Layer  1 

• 

TOE  FAILURE: 

Reinforcing  Layers 

Factor  of  Safety 

Force  of  Layer 

None 

0.373 

- 

• 

One 

(7'  below  surface) 

0.598 

5000  lbs 

Two 

(4'  and  7'  below  surface) 

0.726 

5000  (each) 

• 

Three  0.748 

(2',  4',  and  7'  below  surface) 

1666 

5000 

lbs  (Layer  3) 
lbs  (Layer  1 
and  2) 

* 


ii 
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lbs  per  foot  of  width  were  analyzed  to  produce  the  strength  required 
for  a  minimum  factor  of  safety  of  1.0  (See  Appendix  B).  Analyses  with 
two  and  three  layers  of  reinforcement  revealed  that  5,000  lbs  force  was 
required  to  achieve  a  factor  of  safety  of  1.0  (See  Appendix  B).  It 
should  be  noted  that  with  three  layers  of  reinforcement,  the  top  layer 
was  modelled  as  mobilizing  one  third  the  strength  due  to  the  minimal 
overburden  pressure.  Also,  there  was  minimal  gain  in  factor  of  safety 
from  the  addition  of  the  third  layer.  The  analyses  of  internal 
stability  revealed  the  layers  of  reinforcement  strength  required  for  a 
factor  of  safety  of  1.0  in  the  external  analyses  were  not  strong  enough 
to  achieve  the  required  minimum  stability  requirements  (See  Appendix  B). 
These  results  indicate  internal  stability  is  more  critical,  however, 
the  risk  of  failure  for  an  embankment  overlying  soft  material  is  more 
likely  to  occur  in  the  soft  subgrade  soil  (3). 
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CHAPTER  III 

NUMERICAL  ANALYSES 

? Two  computer  programs  were  utilized  to  analyze  the  perplexing 
roadway  conditions  of  the  Farmer's  Loop  Road  in  Fairbanks,  Alaska. 
FEADAM84:  A  Computer  Program  for  Finite  Element  Analysis  of  Dams  was 
initially  used  to  model  existing  roadway  embankment  conditions. 

SSTIPN:  Soil  Structure  Interactive  Program  was  then  used  to  model 
conditions  for  a  reinforced  roadway  embankment  constructed  on 
permafrost.^ 

3.1  FEADAM84  Analysis 

FEADAM84  is  an  incremental  finite  element  program  for  two- 
dimensional,  plane  strain  analysis  of  earth  and  rockfill  dams  and 
slopes.  It  calculates  the  stresses,  strains,  and  displacements  due  to 
incremental  embankment  construction  and/or  load  application.  The  non¬ 
linear  and  stress  history  dependent  stress-strain  and  volumetric  strain 
properties  of  soils  are  approximated  using  a  hyperbolic  model  modified 
by  Seed  and  Duncan  ( 1 984 ) .  ....(8) 

This  program  is  a  modified  version  of  the  program  FEADAM  develop¬ 
ed  by  Duncan,  Wong,  and  Ozawa  (1980).  The  original  incremental 
analysis  procedure  was  coded  in  the  program  LSBILD  by  Kulhawy,  Duncan, 
and  Seed  (1969).  The  subsequent  program  ISBILD,  by  Ozawa  and  Duncan 
(1973)  incorporated  (a)  isoparametric  elements  with  incompatible  dis¬ 
placement  modes  developed  by  Wilson  et  al.  (1971),  (b)  a  more  accurate 
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procedure  for  assigning  initial  stresses  to  elements,  and  (c)  more 
efficient  computational  techniques,  including  a  new  equation  solver, 
developed  by  Wilson  (1971).  The  program  FEADAM,  by  Duncan,  Wong,  and 
Ozawa  (1980)  incorporated  (a)  a  new  model  for  stress-dependent  bulk 
moduli,  and  (b)  new  criteria  for  differentiation  between  primary 
loading  behavior  and  stiffer,  elastic  unloading/reloading  behavior. 

The  unloading-reloading  model  was  subsequently  deleted  by  Wong  and  Seed 
(1982),  who  also  modified  the  bulk  modulus  model.  ....(8) 

In  its  present  form,  FEADAM84  incorporates  the  following 
modifications  to  the  earlier  program  FEADAM: 

(a)  A  modified  version  of  the  model  for  unloading/reloading 
behavior  has  been  reinstated  on  an  optional  basis, 

(b)  New  criteria  are  employed  for  determination  of  whether  a 
given  soil  element  is  in  a  state  corresponding  to  primary 
loading  or  unloading-reloading  behavior, 

(c)  A  modified  lower  bound  constraint  is  incorporated  in  the 
stress-dependent  bulk  modulus  model, 

(d)  The  new  program  has  the  capability  to  model  weightless, 
linear  elastic  soil  elements  with  zero  initial  stresses. 

3.1.1  Non-linear  Incremental  Finite  Element  Methodology 

A  successive-increment  procedure  is  used  for  approximating  non¬ 
linear  stress  and  stress  history  dependent  behavior  of  soil,  in  which 
progressive  loading  is  divided  into  a  number  of  small  increments,  and 
the  soil  behavior  is  assumed  to  be  linear  within  each  increment.  The 
modulus  values  used  to  model  each  soil  element  are  re-evaluated  each 
increment  in  accordance  with  (a)  the  stresses  In  thp  element-. ,  and  (b) 
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the  previous  stress  history  of  the  element. 

The  incremental  stress-strain  relationship  for  an  isotropic 
material  under  plane  strain  conditions  may  be  expressed  as  follows: 


A  0 

X 
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where  Aox  and  AOy  are  the  incremental  stresses  in  the  x  and  y 
directions  respectively  and  Axxy  is  the  incremental  shear  stress.  Aex 
and  Aey  are  the  incremental  strains  in  the  x  and  y  directions 
respectively  and  AYXy  is  the  incremental  shear  strain.  E  and  B  are 
Young's  modulus  and  bulk  modulus  respectively.  (1) 

The  non-linear,  stress-dependent  and  stress-history-dependent 
soil  behavior  model  is  a  modified  version  of  the  hyperbolic  stress- 
strain,  strength,  and  volumetric  strain  model  proposed  by  Duncan  et  al. 
(1980).  The  original  (1980)  model  has  been  modified  in  order  to  :  a) 
provide  improved  modelling  of  soil  behavior  during  unloading/reloading, 
b)  eliminate  a  source  of  potential  computational  instability  for  some 
types  of  incremental  loading  paths,  and  c)  provide  improved  modelling 

of  bulk  moduli  at  low  stress  levels  and  high  confining  stress . (8) 

The  original  (1980)  model  assumes  that  stress-strain  curves  for 
soils  at  a  given  confining  stress  (03)  can  be  approximated  as  hyper¬ 
bolas  shown  in  Figure  3.1(a).  The  hyperbola  in  this  figure  can  be 
represented  by  an  equation  of  the  form 
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where  Ej.  is  the  initial  tangent  modulus,  and  (ai  -  03'  ;  is  the 

asymptotic  value  of  deviatoric  stress.  In  order  that  the  hyperbolic 
curve  should  fit  the  actual  stress-strain  curve  of  a  given  soil  as 
closely  as  possible  (01  -  03),^  is  always  greater  than  the  compressive 
strength  of  the  soil  (01  -  03) f.  The  two  values  are  related  by  a 
constant  failure  ration  (Rf)  as 

(01  -  03)f  =  Rf(ai  -  03)^  (5) 


and  Rf  is  typically  between  0.6  to  0.9  for  most  soils.  When  the 
deviatoric  stress  exceeds  (a-;  -  03)^  a  failure  condition  is  assumed  as 
shown  in  Figure  3.1(a).  The  deviatoric  stress  at  failure  is  determined 
by  the  familiar  Mohr-Coulomb  failure  criteria  as 


(o,  -  «3>f  - 


cos<)>  +  2o^'  sin$ 
1  -  sin<{> 


where  c  and  <|>  are  the  cohesion  intercept  and  friction  angle  of  the  soil. 
The  model  allows  variation  of  $  as  a  function  of  0 3  as 

°3 

<p  =  <j>o  -  a<J>  •  logio(p  )  (7) 

a 

where  <p0  =  the  soil  friction  angle  at  confining  stress  of  03  =  Pa, 

A$  =  the  reduction  in  <j>  for  each  10-fold  increase  in  03,  and 
Pa  =  atmospheric  pressure,  introduced  in  order  to  make  the  model 
independent  of  the  system  of  units  chosen. 
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The  instantaneous  slope  of  the  hyperbolic  stress-strain  curve 
shown  in  Figure  3.1(a)  is  the  tangent  modulus  Et,  whicn  is  a  function 
of  confining  stress  (03)  and  stress  level  (SL),  and  which  can  be 
expressed  as 


Et  -  Cl  ”  (Rf  *  SL)]2  •  K  •  Pa  •  (o3/Pa)n  (8) 

where  K,n  =  Model  parameters  (constants)  relating  to  the  initial 
modulus  (Ei,  see  Figure  3.1(a))  to  the  confining  stress  o3  as 

Ei  •  K  •  Pa  •  (o3/Pa)n  (9) 

and  SL  =*  Stress  Level,  defined  as  the  ratio 

SL  =  (oi  -  03)/(oi  -  o3)f  (10) 

The  tangent  modulus  modelled  according  to  Equation  8  increases 
with  increasing  confining  stress  (03)  and  decreases  with  increasing 
stress  level  (SL).  ....(8) 

When  the  current  stress  level  was  less  than  the  previous  maximum 
stress  level  (SLmax  past)*  the  unmodified  (1980)  model  assumed  the  soil 
to  be  no  longer  in  a  state  of  primary  loading,  but  rather  in  an 
unloading-reloading  state.  Unloading  and  reloading  were  modelled  as 
linear  and  elastic  as  shown  in  Figure  3.1(b).  The  unloading-reloading 
modulus  were  modelled  as  a  function  only  of  03  according  to  the 
equation 
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Eur  ■  Kur  *  Pa(o3/Pa)n  (11) 

where  Kur  is  typically  1.2  to  3  times  greater  than  K  (the  modulus 
parameter  determining  E^).  ....(8) 

The  bulk  modulus  of  the  soil  is  assumed  to  be  independent  of 
stress  level,  and  can  be  expressed  as  a  function  only  of  03  as 

B  =  KB  •  Pa  *  (a3/Pa)m  (12) 

where  Kg  and  m  are  dimensionless  parameters  (constants).  Modelling  the 
bulk  modulus  (B)  as  being  independent  of  stress  level  results  in 
modelling  an  increase  in  Poisson's  ratio  (v)  with  increasing  stress 
level  because  the  soil  modulus  (E^)  decreases  with  increasing  stress 
level,  and  Poisson's  ratio  can  be  expressed  as 

v  =  1/2  -  Et/6B  (13) 

The  bulk  modulus  in  the  (1980)  model  was  at  all  times  constrained  such 
that  Et/3  <  B  >  17  E^.,  which  in  effect  constrained  Poisson's  ratio  such 
chat  0.0  <  v  >  0.49  (see  Equation  13).  ....(8) 

Together,  the  soil  modulus  (Ej.  or  Eur)  and  the  bulk  modulus  (B) 
define  the  stress-strain  and  volumetric  strain  behavior  of  a  given 
material.  This  relatively  simple,  straightforward  hyperbolic  model  may 
be  formulated  and  applied  in  terms  of  either  total  stress  (0)  or 
effective  stress  (o')  by  using  input  parameters  appropriate  either  for 
total  or  effective  stress  soil  behavior. 


The  unmodified  (1980)  model  suffered  from  two  problems  concerning 
unloading-reload i ng . 

The  first  of  the  two  problems  with  the  previous  model  for  un¬ 
loading-reloading  moduli  is  that  the  simple  stress  level  criteria  used 
for  differentiating  between  primary  loading  and  unloading-reloading 
behavior  does  not  appear  adequate  for  many  purposes.  The  line  of 
constant  stress  level  in  Figure  3.2(a)  represents  the  division, 
according  to  the  earlier  (1980)  model,  between  primary  loading  and 
unloading-reloading  stress  paths  for  an  element  of  soil  currently 
existing  at  a  stress  condition  represented  by  point  A,  where  the  stress 
level  at  A  equals  the  maximum  stress  level  achieved  so  far  (SLmax  pasfc). 
This  division  between  primary  loading  and  unloading-reloading  is 
defined  by  assuming  unloading-reloading  moduli  for  all  conditions  where 
SL  <  Lmax  past*  Subsequent  investigations  have  indicated  that  this 
simple  stress  level  criteria  for  assignment  of  unloading-reloading 
moduli  should  be  modified  in  order  to  include  consideration  of 
variations  in  the  confining  stress  (03)  as  well  as  variations  in  stress 
level.  The  new  (modified)  model  therefore  distinguishes  between 
primary  loading  and  unloading-reloading  based  on  "stress  state"  (SS) 
defined  as 

SS  =  SL  VToTTpj  (14) 

J  o. 

where  SL  =  Stress  level  =  (ai  -  a 3 ) / ( a-|  -  03 ) ,  and 

Pa  =  atmospheric  pressure,  introduced  in  order  to  allow  taking 
the  fourth  power  radical  of  a  dimensionless  number. 
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The  new  model  assumes  primary  loading  moduli  (Et)  are  appropriate  when 
the  current  stress  state  is  greater  than  or  equal  to  the  maximum 
previous  stress  state  (SS  >  SSmax  past)*  and  otherwise  assumes 
unloading-reloading  moduli  (Eur).  Based  on  this  criteria,  the  critical 
stress  level  (SLcrit)  above  which  primary  loading  behavior  is  assumed 
for  a  given  current  confining  stress  (03)  may  be  calculated  as 


SS 


SL 


crit 


max  past 

V^7 ra 


(15) 


Figure  3.2(b)  illustrates  the  division  between  primary  loading 
and  unloading-reloading  stress  paths  based  on  the  new  stress  state 
criteria  for  the  same  initial  conditions  as  were  shown  in  Figure  3.2(a). 
The  stress  space  is  now  divided  by  a  line  which  is  slightly  concave  down¬ 
wards,  and  inclined  slightly  away  from  the  constant  stress  level  line. 

The  second,  and  most  3erious,  problem  with  the  previous  model 
for  unloading-reloading  behavior  is  that  it  can  lead  to  computational 
instability  during  incremental  analyses.  Figure  3.3(a)  illustrates  the 
soil  moduli  resulting  from  application  of  the  unmodified  1980  model. 

At  all  stress  levels  greater  than  or  equal  to  the  maximum  previous 
stress  level  (SLmax  past)  the  primary  loading  modulus  (Et)  is  applied, 
and  at  all  lower  stress  levels  the  unloading-reloading  modulus  (E^p)  is 
applies.  The  sharp  discontinuity  shown  to  occur  in  this  figure  at 
SLmax  past  can  result  in  a  sudden  change  in  modulus  by  a  factor  of  more 
than  10  or  20  at  thi3  point.  This  can  lead  to  computational  instability 
during  incremental  finite  element  analyses,  as  soil  elements  incremen¬ 
tally  subjected  to  loading  at  fairly  constant  stress  levels  may,  as  a 


result  of  minor  fluctuations  in  stress  level,  cross  back  and  forth 
across  the  discontinuity  experiencing  large  and  unwarranted  sudden 
changes  in  modulus.  Attempts  to  minimize  the  impact  of  this  source  of 
potential  computational  instability  by  setting  the  boundary  between 
primary  loading  and  unloading-reloading  at  0.95  SLmax  pa3fc  were 
unsuccessful.  ....(8) 

In  order  to  prevent  this  potential  source  of  computational  in¬ 
stability,  the  new  model  provides  for  a  less  abrupt  transition  between 
unloading-reloading  and  primary  loading  moduli  as  illustrated  in  Figure 
3.3(b),  and  described  below: 

1 .  The  transition  from  unloading-reloading  to  primary  loading 
is  based  on  the  new  stress  state  criteria  as  described 
previously.  Knowing  the  stress  state  (SSmax  past)  which 
defines  the  boundary  between  primary  loading  and  unloading¬ 
reloading  behavior,  the  stress  level  (SLcrit)  at  which 
primary  loading  begins  for  a  given  confining  stress  (03)  can 
be  determined  using  Eq.  15. 

2.  When  SL  >  SLcrit,  the  primary  loading  modulus  (Et)  as 
determined  by  Eq.  8  is  used. 

3.  When  SL  <  3/4  SLcrifc,  the  unloading-reloading  modulus  (Eur) 
as  determined  by  Eq.  11  is  used. 

4.  When  3/4  SLCPit  <  SL  <  SLCPit,  the  modulus  used  is  derived 
by  linear  interpolation  between  EUP  at  3/4  SLCPit  and  Et  at 
SLcrit>  as  shown  in  Figure  3.3(b). 

This  continuous  linear  transition  between  Eur  and  Et  eliminates 
the  sudden  Et/EUP  discontinuity,  and  thus  reduces  the  associated 
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potential  for  computational  instability  inherent  in  the  previous  model. 
The  price  of  this  elimination  of  a  source  of  computational  instability 
is  that  the  new  model  underestimates,  to  some  extent,  the  soil  modulus 
during  the  initial  stages  of  the  transition  from  primary  loading  to 
unloading  behavior.  This  underestimation  is,  however,  of  much  less 
significance  than  the  previous  instability.  In  addition,  the  new  model 
(a)  provides  (correctly)  for  a  modulus  in  the  range  3/^  SLcrit  <  SL  < 
SL-crit  which  is  significantly  higher  than  Et  during  both  unloading  and 
reloading,  (b)  provides  an  appropriate  modulus  (Eur)  in  the  range  SL  < 
3/14  SLcrit,  (c)  provides  an  appropriate  modulus  (Efc)  in  the  range  SL  > 
SLcrit»  and  (d)  provides  an  improved  model  for  soil  behavior  during  the 
transition  from  reloading  to  primary  loading.  ....(8) 

The  previous  (1980)  model  had  a  tendency,  under  certain 
conditions,  to  underestimate  the  bulk  modulus  (B)  and  thus  to  under¬ 
estimate  minor  principal  stresses  as  well  as  soil  strength  and  stiff¬ 
ness.  This  problem  occurred  most  frequently  in  elements  with  very  low 
stress  levels  and  high  confining  stresses.  In  order  to  overcome 
tendency  of  the  bulk  modulus  model  to  underestimate  Poisson's  ratio  at 
low  stress  levels  for  some  soils,  the  lower  bound  constraint  on  the  new 
model  has  been  modified  such  that 

B  .  >  (E./3)  |  2  ~:'S^  }  for  4.  >  2.3°  (16) 

mm  t  1  sin<j>  1  Y 

and  Bmin  =  17  Et  for  <j>  <  2.3°  (16a) 
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This  supercedes  the  unmodified  (1980)  lower  bound  constraint  of 
Bmin  >  Et/3  which  constrained  only  against  negative  values  of  Poisson's 
ratio.  When  <j>  <  2.3°,  is  thus  set  equal  to  0.49,  as  is  appropriate 
for  a  cohesive  soil  under  undrained  conditions.  For  $  >  2.3°,  this 
simple  lower  bound  constraint  on  B  has  the  effect  of  constraining 
Poisson's  ratio  be  greater  than  or  equal  to  vmin  where 


1  ~  s  inct> 
min  2  -  sin<j> 


(17) 


3.1.2  Analysis  Procedure 

The  program  FEADAM84  calculates  the  stresses,  strains,  and  dis¬ 
placements  in  embankments,  simulating  the  actual  sequence  of 
construction  operations  and  post-construction  loading  conditions.  The 
non-linear  and  stress  dependent  stress-strain  properties  of  the  soils 
are  approximated  by  performing  the  analysis  in  increments. 

An  increment  may  consist  of  the  placement  of  a  new  layer  on  the 
embankment,  or  of  application  of  loads  to  a  complete  embankment.  Each 
increment  is  analyzed  twice  (a  two-iteration  process),  the  first  time 
using  modulus  values  based  on  the  stresses  at  the  beginning  of  the 
increment,  and  the  second  time  using  modulus  values  based  on  the 
average  stresses  during  the  increment.  The  changes  in  stress,  strain, 
and  displacement  calculated  during  the  second  iteration  of  each  incre¬ 
ment  are  added  to  the  stresses,  strains,  and  displacements  at  the 
beginning  of  the  increment  to  give  the  total  or  cumulative  values  up  to 


that  stage  of  the  analysis. 


....(8) 


During  each  increment,  each  element  is  checked  to  determine  if 
it  is  in  a  state  of  primary  loading,  unloading/reloading,  tensile 
failure,  or  shear  failure,  and  is  modelled  accordingly  as  follows: 

a)  Primary  Loading.  Primary  loading  moduli  are  used  when  the 

stress  state  of  an  element  is  greater  than  the  maximum 
previous  stress  state  of  the  element,  and  the  minor 
principal  stress  (03)  is  positive.  If  03  is  positive  but 
less  than  0.05  times  atmospheric  pressure,  the  modulus 
values  are  computed  using  03  equal  to  0.05  atmospheric 
pressure.  If  no  unloading-reloading  elastic  modulus  number 
is  specified  for  a  given  soil  type,  then  all  elements  of 
that  type  are  modelled  with  primary  loading  moduli 
regardless  of  stress  state,  unless  the  element  is  in  a  state 
of  tensile  or  shear  failure.  ....(8) 

b)  Elastic  Unloading.  Unloading-reloading  modulus  values  (Eur) 

are  used  when  the  current  stress  state  of  an  element  falls 
below  the  previous  maximum  stress  state.  Linear  elastic 
(stress-state-independent)  materials  may  also  be  modelled 
using  the  program  FEADAM84,  and  these  are  never  considered 
to  be  in  an  unloading-reloading  condition.  ....(8) 

c)  Tensile  Failure.  Tensile  failure  occurs  when  0 3  is  negative. 
During  the  first  iteration,  the  assigned  bulk  modulus  (B)  is 
computed  using  03  equal  to  0.1  times  atmospheric  pressure 
and  the  soil  modulus  (E)  is  set  equal  to  one-tenth  of  the 
bulk  modulus.  If  at  the  end  of  the  first  iteration  the 
element  is  still  undergoing  tensile  failure,  both  of  the 
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moduli  (E  and  B)  are  reduced  to  one- tenth  of  the  values  used 


in  the  first  iteration.  Linear  elastic  materials  are  not 


subject  to  tensile  failure. 


....(8) 


d)  Shear  Failure.  Shear  failure  occurs  when  the  stress  level 


of  an  element  exceeds  95?  of  its  shear  strength.  Modulus 


values  are  modelled  using  Equation  6,  but  corresponding  to  a 


condition  where  the  product  of  stress  level  (SL)  times  the 


failure  ration  (Rf)  equals  0.95,  resulting  in  very  low 


elastic  moduli.  Linear  elastic  materials  are  not  subject  to 


shear  failure. 


....(8) 


3.1.3  Initial  Analysis 


FEADAM84,  in  its  present  form,  is  not  capable  of  calculating 


internal  forces  and  displacements  of  structural  elements  such  as  a 


geotextile  placed  in  a  roadway  embankment.  However,  FEADAM84  is 


capable  of  modelling  the  failure  mechanisms  and  stress  distributions  of 


existing  roadway  embankments  constructed  on  permafrost.  Analysis  of 


geotextile  reinforced  roadway  embankments  will  be  accomplished  with  the 


use  of  SSTIPN  obtained  from  Dr.  J.  M.  Duncan  at  Virginia  Tech.  (See 


Appendix  C) 


3. 1.3.1  Summary  of  Analysis.  Two  material  properties  can  be 


specified  within  the  FEADAM84  program:  (1)  a  linear  elastic  material 


which  is  not  subject  to  shear  failure  and  which  can  take  tension,  and 


(2)  a  stress  dependent  non-linear  elasto-plastic  material  which  obeys 


the  hyberpolic  stress-strain  and  volumetric  relationships  described  in 


Section  3.1 .1 
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Linear-elastic  materials  are  not  applicable  to  this  analysis  and 
will  not  be  discussed.  For  the  stress  dependent,  non-linear  elasto- 
plastic  soil  materials,  a  total  of  nine  parameters  are  required  as 
summarized  in  Table  3.1. 


Three  steps  are  required  to  accurately  model  the  existing 
failure  mode  of  the  Farmer's  Loop  Road  in  Fairbanks,  Alaska.  First, 
the  question  of  "What  were  the  causes  of  roadway  embankment  failure?" 
must  be  answered.  Second,  an  appropriate  representation  of  the  roadway 
embankment  geometry  including  foundation  soil  must  be  defined.  Third, 
an  accurate  representation  of  the  roadway  embankment  performance 
through  the  selection  of  soil  parameter  values  to  be  input  into  the 
program. 

The  failure  mechanisms  associated  with  roadway  embankments 
overlying  weak  foundations  have  been  described  in  Sections  1.1  and  2.2. 
The  failure  of  the  Farmer's  Loop  Road  is  directly  related  to  the 
seasonal  thawing  of  the  frozen  foundation  soil  underlying  the  side 
slopes  of  the  embankment.  This  degradation  can  be  further  amplified  by 
the  insulating  effect  of  snow  over  the  sideslopes  which  prevents  the 
complete  freeze-back  of  the  active  layer  (13).  The  consequential 
consolidation  of  the  thaw  unstable  permafrost  results  in  differential 
settlement  of  the  roadway  embankment  showing  up  as  longitudinal  cracks 
along  the  wearing  surface  (13,  14),  as  shown  in  Figure  1.1. 

The  roadway  embankment  geometry  selected  for  the  numerical 
analyses  (See  Figure  3.4)  is  taken  from  the  State  of  Alaska  Department 
of  Transportation  and  Public  Facilities  (DOT&PF)  Typical  Cross  Section 
of  Farmer’s  Loop  Road  (See  Figure  3.5)  and  the  Engineering  Geology  and 
Soils  Report  (1985).  The  actual  roadway  embankment  geometry  varies 
along  the  length  of  the  road  ana  is  a  function  of  the  amount  of  cut  and 
fill,  grade,  and  superelevation  of  curves.  The  geometry  selected  in 
the  analysis  represents  a  typical  section  along  the  road. 
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Modelling  embankment  performance  is  primarily  dependent  upon  the 
proper  selection  of  soil  parameter  values  to  be  evaluated  by  the 
program.  The  soil  parameters  used  in  the  hyperbolic  model  may  be 
obtained  from  triaxial  compression  tests.  Three  soil  types  are  used  in 
modelling  embankment  performance  in  both  FEADAM84  and  SSTIPN.  The 
roadway  embankment  is  modelled  a3  one  type  for  convenient  use  in  the 
program.  Permafrost  and  thaw-unstable  permafrost  parametric  values  are 
estimated  from  Duncan  et  al.  (1974)  due  to  insufficient  availability  of 
test  data  on  these  types  of  soil  materials.  Since  actual  data  are  not 
available  for  the  thaw-ustable  permafrost,  the  objective  is  to  obtain 
relative  strength  data  for  the  three  soil  types  used  to  model  existing 
conditions.  The  three  soil  types  modelled  are  crushed  rock  for  the 
embankment,  soft  clay  for  the  thaw  unstable  permafrost,  and  hard  clay 
for  the  permafrost.  The  parametric  values  obtained  for  the  crushed 
rock  were  obtained  from  Schauz  (1981).  The  data  obtained  for  soft  and 
hard  clay  were  obtained  from  Duncan  et  al.  (1974).  The  values  obtained 
for  all  three  soil  types  are  summarized  in  Table  3.2. 


Table  3.2  Soil  Parametric  Values  Used  in  FEADAM84  and  SSTIPN 


Parameter 

Crushed  Rock 

Soft  Clay 

Hard  Clay 

K 

6000.0 

40.0 

1000.0 

Kur 

8000.0 

80.0 

1500.0 

n 

0.5 

0.3 

0.4 

Rf 

0.5 

0.9 

0.7 

Kb 

1500.0 

20.0 

500.0 

m 

0.6 

0.2 

0.5 

c 

0  psf 

100  psf 

1000  psf 

4> 

35° 

0° 

40° 

A<j> 

0° 

0° 

0° 

Y 

135  pcf 

106  pcf 

120  pcf 

Kq 

0.5 

0.5 

0.5 

The  relative  locations  of  the  soil  types  and  mesh  used  in  the 
finite  element  analysis  are  shown  in  Figure  3.6.  Location  of  the  talik 
zone  used  in  the  analyses  is  consistent  with  data  obtained  by  the  State 
of  Alaska  DOT&PF  (13,  14,  19).  For  most  cases  reasonably  accurate 
results  can  be  achieved  by  using  eight  or  more  layers  in  the  mesh,  and 
adding  one  layer  at  a  time  to  the  embankment.  However,  six  layers  were 
used  to  investigate  the  relative  accuracy  of  the  model  with  that  of 
Figure  1.1.  Six  layers  also  allowed  more  flexibility  in  areas  which 
required  some  refinement  (See  Figure  3.7).  Figure  3.8  represents  the 
deformed  state  of  the  original  model.  The  failure  mode  associated  with 
the  Farmer's  Loop  Road  can  be  detected  in  Figure  3.8.  The  deformed 
geometry  plot  is  performed  using  Program  GRAPH  (See  Appendix  D).  The 
consolidation  of  talik  zone  and  severe  surface  deformation  of  the  top 
of  the  embankment  are  clear.  Through  analysis  of  this  figure,  the 
material  properties  are  reasonably  accurate  for  further  use  in  SST1PN. 

3.?  SSTIPN  Analysis 

SST1PN  calculates  stresses,  strain  and  displacements  in  soil 
elements  and  internal  forces  and  displacements  in  structural  elements 
by  means  of  analyses  which  simulate  the  actual  sequence  of  construction. 
The  non-linear  and  stress-dependent  stress-strain  properties  of  the 
soils  are  approximated  using  the  procedures  developed  by  Kulhawy, 

Duncan,  and  Seed  (1969).  The  structural  materials  are  assumed  to 
behave  linearly.  The  soil  model  used  is  predecessor  to  the  existing 
soil  model  used  in  FEADAM84.  The  purpose  of  utilizing  SSTIPN  is  to 
analyze  the  effects  of  placing  different  structural  materials  within  a 


FEADAM84  EMBANKMENT  MODEL 


roadway  embankment  constructed  on  permafrost. 

3.2.1  Assumptions 

The  results  of  FEADAM84  reveal  some  important  points.  First, 
the  soil  parameters  selected  for  the  analyses  are  appropriate  and 
accurately  model  roadway  embankment  performance  in  permafrost  regions. 
The  deformed  roadway  embankment  model,  Figure  3*8,  depicts  the 
distressed  condition  as  a  result  of  thaw-unstable  foundation  soil 
located  beneath  the  side  slopes  of  the  embankment.  Second,  the 
geometry  and  location  of  the  soil  types  are  relatively  accurate  within 
the  model.  The  relative  location  of  the  soil  material  properties, 
Figure  3 - 6 ,  allows  the  roadway  embankment  model  to  behave  consistently 
with  DOT  &  PF  articles  (9,  21). 

The  conditions  of  the  SSTIPN  analyses  are  as  follows: 

1 .  Embankment  geometry  is  identical  to  both  STAB CM  and  FEADAM84 
analyses.  Foundation  soil  height  =  10  feet.  Embankment 
soil  height  *  7  feet.  Total  number  of  soil  elements  if 
increased  from  48  to  78  elements. 

2.  Five  different  strength  modulae  are  analyzed  for  the 
geotextile  properties.  In  addition,  five  different 
conditions  are  analyzed  for  each  separate  geotextile 
strength  modulus. 

3.  Soil/geotextile  interface  is  modelled  for  no  slip  between 
soil  elements  and  structural  elements. 

Figure  3.9  depicts  the  mesh  that  is  used  in  the  SSTIPN  analyses. 
The  foundation  and  embankment  numbers  of  layers  are  increased  to  three 


and  six,  respectively  to  provide  more  accurate  results  of  displacements 
and  stresses.  Figure  3*10  depicts  the  location  of  the  three  soil 
material  types  in  the  analyses.  The  location  of  the  thaw-unstable 
material  beneath  the  side  slopes  is  consistent  with  articles  from  the 
state  of  Alaska  DOT&PF  (9,  15,  21). 

The  five  geotextile  strength  modulus  are  displayed  in  Figure 
3.11.  Tensile  strengths  of  1000,  2000,  and  4000  pounds  per  inch  with 
strains  of  5  and  20  percent  are  evaluated.  SSTIPN  does  not 
differentiate  between  modulae  with  different  strain.  SSTIPN  calculates 
the  stiffness  of  each  member  in  force  per  unit  width  without  considera¬ 
tion  to  percent  elongation.  Therefore,  a  material  with  strength  of 
1000  pounds  per  inch  at  5  percent  strain  is  identical  to  a  material 
with  a  strength  of  4000  pounds  per  inch  at  20  percent  strain. 

The  five  conditions  analyzed  for  each  strength  modulae  are  as 
follows: 

1.  Three  layers  of  geotextile  are  located  within  the  embankment. 
The  first  layer  located  at  the  embankment/ foundation  inter¬ 
face  and  the  second  and  third  layers  located  three  and  five 
feet  above  the  embankment/ foundation  interface, 
respectively. 

2.  Three  layers  of  geotextile  with  identical  locations  as  the 
first  condition  with  a  pre-tension  of  1000  pounds  per  foot 
on  all  three  layers. 

3.  One  layer  of  geotextile  located  at  the  embankment/foundation 
interface. 


4.  One  layer  of  geotextile  located  at  the  embankment/ foundation 


interface  with  a  pre-tension  of  1000  pounds  per  foot. 

5.  One  layer  of  geotextile  located  at  the  embankment/ foundation 
interface  with  a  pre-tension  of  1000  pounds  per  foot.  The 
geotextile  extends  10  feet  beyond  the  toe  of  the  embankment. 
SSTIPN  has  the  capability  of  modelling  friction  between  soil 
and  structural  elements.  The  geotextile/soil  interface  is 
modelled  with  total  friction,  i.e.,  no  slip  between  the  soil 
and  geotextile.  The  overall  accuracy  increases  only  by 
about  5  percent  with  the  addition  of  friction  elements. 
Modelling  without  friction  elements  is  convenient  for 
repetitive  analyses  and  does  not  reduce  the  accuracy  enough 
to  cause  a  significant  difference. 

3.2.2  Results  of  SSTIPN  Analysis 

The  analysis  of  SSTIPN  results  includes  obtaining  displacement 
and  stress  data  from  the  output  to  evaluate  the  effect  of  each 
different  modulae  given  the  different  conditions.  The  analysis  of  the 
roadway  embankment  with  three  layers  of  geotextile  did  not  provide 
reliable  results  due  to  instability  of  SSTPIN  in  calculating 
displacements.  Therefore,  analyses  of  the  five  modulae  with  the  first 
two  conditions  mentioned  will  no  longer  be  discussed. 

Displacement  data  extracted  from  SSTIPN  output  for  all  five 
different  strength  modulae  with  the  last  three  conditions  are  plotted 
using  the  GRAPH  in  Appendix  D.  Figures  3.12  through  3.27  depict  the 
deformed  geometry  for  each  case.  Figures  3.28  through  3.^3  the 


deviatoric  stress  distribution  for  the  identical  conditions  used  to 


plot  the  deformed  geometry. 
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Conditions:  1  Layer  of  Geotextilc.  Extended  10  Ft 
Seycnd  Toe  of  Emoankment 
Pre-tension  =  1000  lbs/ft 
Geotextile  Strength:  1000  lbs/ in  @  t  =  20 % 
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Conditions:  1  Layer  of  Geotextile  Extended  10  Ft 


FIG.  3.21  DEFORMED  EMBANKMENT  GEOMETRY  OF  SSTIPN  ANALYSIS 


'RMED  EMBANKMENT  GEOMETRY  OF  SSTIPN  ANALYSIS 
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FIG.  3.29  DEVIATORIC  STRESS  CONTOUR  OF  SSTIPN  ANALYSIS 
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Conditions:  1  Layer  of  Geotextile 

No  Pre-tension 

Geotextile  Strength:  4000  lbs/in  @ 


3.2.2  Conclusions  and  Recommendations 

There  are  two  major  factors  which  determine  whether  or  not  a 
numerical  analysis  is  reliable.  First,  existing  conditions  must  be 
modelled  properly.  Geometry  and  material  properties  must  be  as  close 
to  actual  conditions  as  possible  to  accurately  simulate  behavior  of  the 
model.  Second,  the  program  must  be  able  to  evaluate  conditions  which 
may  exist  in  the  field.  The  soil  model  must  be  able  to  determine 
appropriate  values  for  strength  to  simulate  actual  conditions.  The 
hyperbolic  model  used  in  SSTIPN  models  soil  conditions  very  accurately 
up  to  the  point  of  failure.  Once  failure  of  the  soil  results,  the 
model  i3  not  capable  of  simulating  accurate  conditions.  Appendix  E 
contains  detail  results  of  all  analyses.  Close  observation  of  the 
analyses  reveals  the  occurrence  of  shear  and  tensile  failure  within 
soil  elements  located  within  the  thaw-unstable  soil  and  roadway  embank¬ 
ment.  Failure  occurs  within  the  foundation  soil  from  the  placement  of 
the  embankment's  first  layer.  With  the  placement  of  the  additional 
five  layers  of  embankment  soil,  an  increasing  number  of  soil  elements 
enter  a  failure  mode.  Observation  of  stress  contours  for  all 
conditions  do  not  reveal  a  major  redistribution  in  stress  throughout 
the  embankment  and  foundation  as  a  result  of  the  addition  of  a 
geotextile  and  increases  of  strength  in  the  geotextile.  The  stress 
contours  of  the  analysis  with  no  layers  of  geotextile  is  not  very 
different  from  the  analysis  with  one  layer  of  geotextile  with  a 
strength  of  4000  pounds  per  inch  at  5  percent  strain.  The  inability  of 
the  hyperbolic  model  to  accurately  simulate  soil  conditions  after 
failure  is  the  primary  cause  for  the  difference.  With  soil  elements  in 
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a  failure  mode  throughout  each  analysis,  the  computation  of  stresses 
within  the  embankment  apparently  do  not  differ  very  much.  In  the  final 
analysis,  the  resulting  stress  distributions  are  similar. 

Observation  of  the  deformed  roadway  embankment  models,  Figures 
3.12  through  3-27,  reveal  subtle  differences  between  each  condition. 

The  inability  of  the  hyperbolic  model  to  simulate  soil  properties  after 
failure  has  the  same  effect  on  displacements  as  it  does  on  stress 
computation.  However,  the  diagrams  display  more  displacement  with 
weaker  geotextile  properties.  This  indicates  the  positive  effect  of 
stronger  geotextile  properties.  It  is  difficult  to  assess  which 
property  has  the  most  positive  effect  due  to  the  programs  inability  to 
distinguish  between  different  strain  values.  However,  geotextiles 
which  have  higher  strength  and  higher  modulae  tend  to  maintain 
embankment  integrity  and  reduce  the  failure  mechanisms  associated  with 
embankment  stability.  The  more  the  embankment  is  allowed  to  deform 
over  the  thaw-unstable  region,  the  more  the  risk  for  longitudinal 
cracking  in  the  wearing  surface. 

The  addition  of  a  soil  berm  and  geotextile  extension  outside  the 
toe  of  the  embankment  causes  the  degradation  of  thawing  to  occur 
beneath  the  berm  more  than  the  sideslopes  of  the  embankment.  The 
result  is  consolidation  of  the  thaw-unstable  soil  beneath  the  berm  more 
than  beneath  the  sideslopes  of  the  embankment.  This  increases 
embankment  stability  and  reduces  the  risk  of  longitudinal  cracking  in 
the  wearing  surface. 


Through  further  observation  of  geotextile  forces  in  Appendix  E, 
the  failure  capacity  is  not  attained  in  any  analysis.  This  indicates 


that  geotextile  integrity  is  not  impaired  and  will  not  fail  under  these 
conditions. 

A  comparison  of  Figure  3.8  and  3.12  reveals  a  discrepancy 
between  the  resulting  displacements  in  the  FEADAM84  and  SSTIPN  analysis. 
The  soil  model  employed  in  the  SSTIPN  program  is  an  older  version  of 
the  model  used  in  the  FEADAM84  program.  Displacements  of  the 
embankment  model  should  have  been  similar  given  the  same  condition. 

The  FEADAM84  deformed  geometry,  Figure  3*8,  shows  more  displacement 
than  the  SSTIPN  deformed  geometry,  Figure  3.12.  The  incr  iased  number 
of  soil  elements  in  the  SSTIPN  analysis  should  add  more  flexibility  to 
the  model.  From  the  results  obtained,  the  opposite  occurs.  The  cause 
of  this  discrepancy  is  primarily  attributed  to  the  different  generation 
of  hyperbolic  soil  model. 

Also,  proper  modelling  of  soil  conditions  after  failure  is  best 
achieved  with  a  plasticity  model.  The  addition  of  this  type  of  soil 
model  in  this  program  would  greatly  enhance  its  effectiveness  in 
evaluating  roadway  embankments  constructed  over  weak  soil. 
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APPENDIX  A 

STABGM  INITIAL  ANALYSIS 

Figure  2.4  depicts  the  conditions  of  the  initial  STABGM  analysis. 
Detailed  results  of  the  initial  analyses  are  contained  herein.  Minimum 
factor  of  safety  data  from  analyses  with  no  layers  of  reinforcement 
with  different  foundation  geometry  are  enclosed.  Comparison  of  minimum 
factor  of  safety  data  with  three  and  four  layers  of  reinforcement  are 
also  enclosed. 
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Program  STABR  —  ':.-rsicn  2.84  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
nerma-frost  depth  5  feet,  10  -Feet  -from  centerline 


Control  data 

NUMBER  OF  SPECIFIED  CENTERS  0 
NUMBER  OF  DEPTH  LIMITING  TANGENTS  3 
NUMBER  OF  VERTICAL  SECTIONS  12 
NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  PORE  PRESSURE  LINES  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  -  c/-*.  40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  MET  IOP 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF  2.0 

ALL  CIRCLES  TANGENT  TO  DEPTH,  2B.0,  30.0,  35.0, 

GEOMETRY 
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.  0 

10.  0 

13.0 

15.0 

20 . 0 

32.  0 

42.0 

64 . 0 
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7  4.0 

T.  CRACKS 

20 . 0 

20 . 0 

20 . 0 

20 . 0 

20 . 0 

20 . 0 

25.0 

25.0 

25 . 0 

25.0 

25.0 

W  IN  CRACK 

20 .  0 

20 .  0 

20.  0 

20.0 
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20 , 0 

25 . 0 

25.0 

25.  :> 

25.  0 

25.0 

BOUNDARY  1 

20 . 0 

20 . 0 

20.0 

20 . 0 

20.0 

20 . 0 

25 . 0 

23 . 0 

25 . 0 

25 , 0 

25.  C 
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23 . 0 

25 . 0 

25.0 

25.0 

25.  0 

25 . 0 

25 . 0 
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25.  0 

25.0 

25 .  0 
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35 . 0 
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REINFORCING  FORCE  DATA  AT  0  LEVEL  1  £ ) 
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BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  5  feet,  10  feet  from  centerline 
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TANGENT 

RADIUS  CO  CENTER 
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FS (BISHOP 
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Program  STABR  —  Version  2. 34  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
oarmaf rest  depth  5  reet.  15  feet  -from  csntierlins 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  0 
NUMBER  OF  DEPTH  LIMITING  TANGENTS  3 
NUMBER  OF  VERTICAL  SECTIONS  12 
NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  PORE  PRESSURE  LINES  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF  2. 

ALL  CIRCLES  TANGENT  TO  DEPTH,  28.0,  30.0,  35.0. 

GEOMETRY 


SECTIONS 

.  0 

15.0 

18.0 

20.0 

25.  0 

32.0 

42 . 0 

59.0 

64 . 0 

66.0 

69.0 

T.  CRAChS 

20 . 0 

20.0 

20.0 

20.0 

20.0 

20.0 

25 , 0 

25.0 

Ci 

23 . 0 

25 . 0 

W  IN  CRACt" 

20 . 0 

20 . 0 

20 . 0 

20.  0 

20 . 0 

20 . 0 

25.  0 

25.0 

25.0 

25.  0 

25 .  0 

BOUNDARY  i 

20.0 

20 . 0 

20. 

20.0 

20,0 

20 . 0 

25.  0 

25.0 

25 . 0 

25 . 0 

25 . 0 

BOUNDARY  2 

25.0 

25.0 

25.  0 

25 . 0 

25.0 

25.  0 

25.0 

25.0 

23.0 

25.  0 

25.0 
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35.0 

35 . 0 
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25.0 

25 . 0 

25.0 

25.0 

25 . 0 

25 , 0 
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35.0 
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no  r, 
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35 . 0 
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35.0 

35 . 0 

35.0 

35 . 0 

35 . 0 
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>'0  «  . 

35 . 0 
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permafrost  depth  5  -feet,  15  -feet  -from  centerline 
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p 
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C7  ^  1 

Q 
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10 

23 . 0 
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33 . 0 
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.  950 
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23 . 0 
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.  950 
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28 . 0 
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MINIMUM" 


966  FOR  THE  CIRCLE  OF  CENTER 


i  8 . 0 ) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
psrmdTrost  depth  5  -feet,  15  -feet  from  centerline 
NUMBER  TANGENT  RADIUS  (X)  CENTER  Cv>  CENTER  FE'BI 
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874  FOR  THE  CIRCLE  OF  CENTER  (  37.0,  18.0) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  5  -feet,  15  -feet  from  centerline 

NUMBER  TANGENT  RADIUS  <X;  CENTER  ’Y)  CENTER  FS  < BISHOP '  P'i  QMS, 
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Program  S I ABR  —  Version  2.84  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  5  -faet,  20  test  -from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


i  a 


0 

0 


SEISMIC  COEFFICIENT  SI, SI 
UNIT  WEIGHT  OF  WATER 


. 00 ,  . 00 
62.40 


SEARCH  IS  BASED  GN  BISHOP  MODIFIED  METHOD 


SEARCH  STARTS  AT  CENTER  <  37.0,  22.0)  WITH  FINAL  GRID  OF  2.0 

ALL  CIRCLES  TANGENT  TO  DEPTH,  28.0,  30.0,  35.0, 

GEOMETRY 


SECTIONS 

.0 

20.  0 

23.0 

25. 0 

30.0 

32.  0 

42.0 

54 . 0 

5? .  0 

6  1 . 0 

64.0 

T.  CRACKS 

20.0 

20 . 0 

20 . 0 

20.0 

20.  0 

20. 0 

25 . 0 

25.0 

25. 0 
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25.0 

W  IN  CRACK 

20.0 

20.0 

20.  0 
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20.  0 

20 . 0 

25 . 0 

25 . 0 

25.0 

25.0 

25.0 
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20 . 0 

20 . 0 
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20 . 0 

20.0 
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25 . 0 

25 . 0 
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25.0 
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i Program  STABR  —  Version 


ini 


-DOS' 


0  BISHOP  MODIFIED 

Opermaf cost  depth  5  feet 


AND/OR  ORDINARY  METHOD  OF 
20  feet  -from  centerline 


SLICES 
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SPECIFIED  CENTERS 
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PORE  PRESSURE  LINES 
POINTS  DEFINING  COHESION  PROFILE 
=  . 00 ,  . 00 
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25.0 

25,0  25.0  05.0  25, 

BOUNDARY  2 

25.0  25.0  25.0  25.0  25.0  25.0 

25 « O 

25.0  25.0  25.0  25. 

BOUNDARY  3 

35.0  35.0  25.0  25 .  C  25..  0  25.''- 

25 ,.  0 

25.0  25-0  25.0  35.. 

BOUNDARY  4 

35 . 0  35 .  < 

<  28 . 0  30 . 0  35 . 0  35 . 0 

35 . 0 

35.0  30. 0  28.0  35. 

BOUNDARY  5 

35,0  35.0  35.0  35.  o  35.0  33,.  O 

35 . 0 

35 . 0  35 . 0  35 . 0  33 „ 

ORE INFORCING  FORCE  DATA  AT 

3  LEVEL (s) 

0  Y~  20.50 

NO.  OF  FORCE  POINTS=  4 

X 

FORCE 

i", 

.  0 

30.  0 

1 666 . 0 

3.0 

1.666,  0 

„  0 

.  0 

0  Y“  22.00 

NO.  OF  FORCE  POINTS^  4 

X 

FORCE 

36 . 0 

,  ( 

) 

33-  0 

5000 .  < 

\ 

3  „  0 

S0O0 .  ( 

) 

i"i 

.  ( 

) 

0  Y-  24.00 

NO.  OF  FORCE  P0IN7S=  4 

X 

FORCE 

40 .  0 

f 

37 .  0 

5000. ' 

) 

3 . 0 
» 0 

5i.i00 . 

.  ( 

\ 

030 IL  PROPERTIES 

0  L.AYER 

COHESION  FRICTION  ANGLE 

DENS.1  T 

/ 

i 

I 

*  ZC')  „ 

■L  -«V  ,\ 

■ 

1 000 ,  < 

1  40.0 

1 2«“> 

-■ 

1  AO ,.  > 

.  O1 

!  05  ,  0 

l 

i 0O0 . 1 

•  20.  o 

0  BISHOP  MODIFIED  AI-!D/Op  ORDINARY  METHOD  OF 

Ti.  rOES 

'.'perma-f ro?t  depth 

5  feet,  20  feet  from  centerline 

ONUMBER  TANGEN  T 

RADIUS 

■'.<)  CEN'EP  ' V  i  CEN-Ei' 

r  r. ;  r  i  s 

HOF;  FS  1  (.'MS *  iFS  ■  BE 

i  2B.0 

6 . 0 

-•7.0  DO.,*!1 

cfim  r:: 

R  BFLOW  T  I'FEPPOLATED  O 

2  28. 0 

6  •  0 

33  o  72,'.' 

LvL'.P'i  I 

F  BEt.uW  JNTEFFOLATED  f. 

l? 

I 


TS  CD 


!T:^LTyyA,lt\r*yvi'T,ir«VTJ.^l':i(l.’WVV  trcgTOmff^H  l~v:.-vi-y.L~>.'-^  im  ir*V».vr>  lTfc.U-tunA7W.tWH.1* 


Program  STABP  —  Version  2.34  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  5  -feat,  20  feet  from  centerline 


h 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


O 


12 

5 

0 

0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37,0,  22.0)  WITH  FINAL  GRID  OF 

ALL  CIRCLES  TANGENT  TO  DEPTH,  28.0,  30.0,  35.0, 

GEOMETRY 


,  0 


SECTIONS 

.  0 

20 . 0 

23.0 

25 . 0 

30 . 0 

32 . 0 

42.0 

54 . 0 

59.0 

61.0 

64.0 

T.  CRACKS 

20,  0 

20.  0 

20.0 

20.0 

20 . 0 

20 . 0 

25 . 0 

25.0 

25.0 

25 . 0 

25.0 

W  IN  CRACK 

20.  0 

20. 0 

20 . 0 

20.  0 

20 . 0 

20 . 0 

25.0 

25 . 0 

25 . 0 

25 . 0 

25,0 

BOUNDARY 

J 

20 . 0 

20.  0 

20 . 0 

20 . 0 

20 . 0 

20 . 0 

25 . 0 

25 . 0 

25.0 

25.0 

25 . 0 

BOUNDARY 

n 

25.0 

25.  0 

25 . 0 

25 . 0 

25.0 

25.0 

25 . 0 

25 . 0 

25.0 

25.0 

25 . 0 

BOUNDARY 

3 

35 . 0 

35 . 0 

25 . 0 

25 . 0 

25.0 

25 . 0 

25 . 0 

25.0 

25 . 0 

25 . 0 

35  „  0 

BOUNDARY 

4 

35 . 0 

35.  0 

28.0 

30 . 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

30.  0 

28.0 

35.  ■') 

BOUNDARY 

5 

35. 0 

{ } 

35.  0 

35.  0 

35.0 

35 . 0 

35 . 

35 . 0 

35.  0 

35 . 0 

35.  0 

REINFORCING  FORCE  DATA  AT  4  LEVEL (s) 


Y  - 


20 . 50 

NG.  OF 

X 

FORCE 

3  3 . 0 

,  0 

30 . 0 

1. 666 . 0 

3. 0 

1 666 . 0 

.  0 

.  0 

22 . 00 

NO.  OF 

x 

FORCF 

36 . 0 

.  0 

33.  0 

3000 . 0 

3 . 0 

5000 . 0 

,  0 

.  'j 

24 . 00 

NG,  OF 

X 

FORCE 

40.  0 

.  0 

37 . 0 

5000. C 

3 . 0 

50 O'"-.  0 

5-  4 


FQFf 


POINTS— 


b 


i.'-  yv  \T*  ’.’V  Vm*  "  V*  .  '  - 


M 


42 . 0 
79  „  0 


5000 .  0 
3000 *  0 

.  0 


L  PROPERTIES 

LAYER 

1 


COHESION 
.  0 

1 000 . 0 
100,0 
1 000 . 0 


FRICTION  ANGLE 
35 . 0 
40.  0 
.  0 
40.  0 


DENSITY 
135.0 
i  20 . 0 
105.  0 
120.0 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
■mafrost  depth  5  -feet,  20  -Feat  from  centerline 
1BER  TANGENT  RADIUS  (X)  CENTER  (Y;  CENTER  FS  f B I 


FSf BISHOP)  FS (QMS)  dFS(BS) 


1 

28.0 

6. 0 

37 . 0 

22 . 0 

28.0 

6.0 

33 . 0 

22 . 0 

PI 

28.0 

1 0 . 0 

37.0 

18.0 

4 

28.0 

6.0 

41.0 

22 , 0 

5 

28.0 

2.0 

37.0 

26 . 0 

6 

28.0 

10.0 

35 . 0 

18.0 

7 

28.  0 

12.0 

37 . 0 

1 6 .  0 

S 

28.0 

10.0 

39.0 

18.0 

o 

/ 

28.  0 

8.0 

37 . 0 

20 . 0 

10 

28 . 0 

8 . 0 

35 . 0 

20 . 0 

u 

28 . 0 

8.  0 

39 . 0 

20 . 0 

12 

28.0 

6.0 

37 . 0 

22 . 0 

13 

28.0 

1 0 . 0 

35.0 

18.0 

14 

28.0 

10.0 

39.0 

18.0 

15 

28.0 

o.  0 

39 . 0 

22 . 0 

16 

28.0 

6 . 0 

35 . 0 

22 . 0 

c 

MINIMUM* 

4.473  FOR 

THE  CIRCLE  OF 

CENTER 

BISHOP 

MODIFIED 

AND/OR  ORDINARY 

'  METHOD 

■rmi 

1  Ivi  l.J  I 

a-frost  depl 

rr-  TAunrrM" 

:h  5  -feet, 

r  DAn  t  i  ic 

20  -feet  from  c 

/  v  \  rtrurrr. 

enter  1 \  i 

CENTER 

CENTER 

4.778 

CENTER 

CIRCLE 

5.081 

4.962 

5.081 

4.473 

4.897 

4.876 

CENTER 

5.081 

5.081 

CENTER 

CENTER 


BELOW  INTERPOLATED  CRN 
BELOW  INTERPOLATED  C r! 

4.67?  3.811 

BELOW  INTERPOLATED  CR£ 
OUTSIDE  SLOPE 


4.965 

4 .  O  4  4 

4 . 383 

3 , 974 

4.965 

4.  045 

4.349 

3.496 

4.722 

3.795 

4.  "27 

3.799 

BELOW  INTERPOLATED  CRf 
4.965  4.044 

4.965  4.045 

BELOW  INTERPOLATED  CR£ 
BELOW  INTERPOLATED  OF; 


1 

30.0 

8.0 

2 

30 .  'I1 

8.0 

pi 

30. 0 

12.0 

4 

30 . 0 

8 . 0 

5 

30 . 0 

4.0 

b 

30 . 0 

1 2 , 0 

7 

70 . 0 

1 4 , 0 

8 

30 .  C- 

I  2 . 0 

9 

30  „  o 

]  0 . 0 

0 

30 0 

1 0 . 0 

i 

30 .  0 

10.0 

TO  a 

8 . 

p 

10  N  0 

12.0 

4 

30  -  0 

1 0 . 0 

5 

30. 0 

8.  0 

6 

30, 0 

!.  2 . 0 

35  -  0 

22 . 0 

CENTER 

BELOW 

ji.U 

22.0 

CENTER 

BELOW 

35 . 0 

18.0 

7  -  328 

39 . 0 

22 ,  o 

CENTER 

BEL  GW 

35 . 0 

26 . 0 

CIRCLE 

ours  I 

33 . 0 

*  f-  „  0 

75  2 

35 . 0 

1 6 . O 

3 . 520 

77T0 

i  8 . 

z, 

35 . 0 

20 . 0 

3.  ■;  i9 

33 , 0 

20 . 0 

3.574 

i't 

'"Vi  i’l 

2.869 

33  -  0 

^  —>  f  , 

CENTER 

Bfli.  OW 

37 . 0 

18.'.' 

-  20o 

39 . 0 

o 

7 .  0  1  9 

37  -  O 

22 .  ’I* 

CE! i'Fr" 

BEL  OvJ 

3,249 


I  C.L.-  Civ 

2. 41S 


3.452  2.606 

2.929  Z. OS 6 

3.,  147HEG.  RESIST. 

~  "oo  i  or; 


1  / 

jO .  (j 

1 2 .  o 

3V .  ' 

18,0 

3.  328 

.  249  2.417 

18 

70 . 0 

8.0 

39.0 

22.0 

CENTER 

BELOW 

INTERPOLATED  CF 

1? 

30 , 0 

8.0 

35,0 

22 . 0 

CENTER 

BELOW 

INTERPOLATED  CF 

F.S. 

MINIMUM- 

2.869  FOR 

THE  CIF<CLE  OF 

CENTER  < 

37 . 0 ,  20 . 

0) 

BISHOP  MODIFIED 

AND/OR  ORDI 

NARY  METHOD  OF 

SLICES 

pormai 

rost  depth 

5  feet , 

20  f  set  f  r 

am  centerline 

NUMBER 

TANGENT 

RAD I  US 

(X)  CENTE 

R  (Y>  CENTER 

FS (BISHOP i  FSCOMS) 

dFS ( B9 ; 

1 

35-  0 

13,0 

35.0 

22.  0 

CENTER 

BELOW  INTERPOLATED  CF 

35.0 

13.0 

31.0 

22 . 0 

CENTER 

BELOW  INTERPOLATED  CP 

•i 

35 . 0 

1 7 . 0 

33 . 0 

18,0 

5.838 

5 . 347 

1 , 008 

4 

35.0 

1 3 . 0 

39.0 

22 . 0 

center 

BELOW  INTERP 

OLATED  CF: 

5 

35.0 

9.0 

35 . 0 

26.0 

CENTER 

BELOW  RIGHT 

INTERSECT 

6 

35  -  0 

17.0 

32 . 0 

18.0 

7 . 489 

6.778 

1 , 027’ 

7 

35  *  0 

19.0 

35 . 0 

i  6 . 0 

6.202 

5.722 

1.148 

8 

35 . 0 

17.0 

37 . 0 

18.0 

2.656 

2.763 

1 . 020 

9 

35 . 0 

15.0 

35,0 

20 . 0 

CT  -r-y-r 

4.837 

795 

1 0 

35 . 0 

19.0 

37.0 

16.0 

4.511 

4 . 30  i 

1 .  1 62 

11 

35 . 0 

17.0 

39.0 

18.0 

1 . 882 

1 , 837 

1 , 03c 

12 

75 . 0 

15.0 

37 . 0 

20 . 0 

1 . 654 

1.64  4 

.  T795 

L3 

35. 0 

15-0 

35 . 0 

20 . 0 

cr  ^ 

4 . 837 

.  795 

14 

75 - 0 

15.0 

39.0 

20 . 0 

1.671 

1.623 

,  6 1 0 

i5 

35.0 

13.0 

37 . 0 

22.  0 

CENTER 

BELOW  INTERPOLATED  CRi 

1  o 

35 . 0 

17.0 

35 . 0 

18.0 

5.838 

5.347 

L .  008 

17 

35.  0 

17.0 

39 . 0 

18.0 

1.882 

1 . 837 

1 . 039 

18 

35.  0 

13.0 

39.  0 

22.  0 

CENTER 

BELOW  INTERPOLATED  CR, 

19 

35 . 0 

13.0 

35 . 0 

22 . 0 

CENTER 

BELOW  INTERF 

■OLATED  CR: 

F.S.  MINIMUM55  1. 

654  FOR 

THE  CIRCLE 

OF  CENTER  ( 

37.0,  20. 

0) 

•  V/ 

$: 

t> 


fei 


APPENDIX  B 


RESULTS  OF  STABGM  ANALYSIS 

Figure  2.5  depicts  the  conditions  of  the  revised  STABGM  analysis. 
Detailed  results  of  the  revised  analysis  are  contained  herein. 

Analysis  with  no  layers  of  reinforcement  display  the  critical  factor  of 
safety  exists  at  a  depth  of  10  feet  below  the  embankment/ foundation 
soil  interface.  Analyses  with  one  layer  of  reinforcement  are  enclosed. 
Strengths  of  5,000,  7,000,  9,000,  and  8,000  pounds  per  foot  were 
analyzed  to  obtain  a  minimum  factor  of  safety  of  1.0  at  a  strength  of 
8,000  pounds  per  foot.  Analyses  with  two  and  three  layers  of 
reinforcement  are  also  enclosed.  Internal  stability  analyses  with 
none,  one,  two,  and  three  layers  of  reinforcement  are  enclosed  for 
comparison  with  external  stability  analyses. 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  0 

NUMBER  OF  DEPTH  L I M I T T NG  TANGENTS  3 

NUMBER  OF  VERTICAL  SECTIONS  12 

NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 

NUMBER  OF  POPE  PRESSURE  LINES  0 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 

SEISMIC  COEFFICIENT  SI ,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  -  62.40 

SEARCH  IS  BASED  O '  BISHOP  MODIFIED  METHOD 

SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  Gi 


ALL  CIRCLES  TANGENT  TO  DEF-'I H ,  28.0,  30.0 


GEOMETRY 

SECTIONS  .0  10.0  13.0  15. 


•35.0 


T.  CRACKS 

18.0 

18.0 

18.0 

W  IN  CRACK 

18.0 

18.0 

18.0 

BOUNDARY  1 

18.0 

5  3.0 

18.0 

BOUNDARY  2 

25.0 

25,0 

Z5.0 

BOUNDARY  3 

35 , 0 

35.0 

25.0 

BOUNDARY  4 

35.  0 

7-5.0 

28.0 

BOUNDARY  5 

35, 0 

35 . 0 

35 . 0 

REINFORCING  FORCE  DATA  AT  0  LEVEL (s) 


SOIL  PROPERTIES 
LAYER 

i 


COHESION  FRICTION  ANCL 
.  0  35 . 0 

1000.0  40..  0 

1  '""I* 1 .  0  .  0 

1000.0  40.0 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHO 
p&rmar^Dst  dc?pth  7  -feet.,  10  ~e-?n  £  •  csm  canter: 


V 

28 .  0 

8.0 

77 . 0 

20. 0 

CENTER  BELOW 

INTERPOLATED  CF.  j& 

10 

28.0 

12.0 

35 , 0 

J  6 . 0 

.  779 

.710  .004 

11 

28.  0 

12.0 

39.0 

16.0 

.  7  79 

.710  .000  H 

12 

28.0 

8.0 

39 . 0 

20 , 0 

CENTER  BELOW 

INTERPOLATED  CF  ££ 
INTERPOLATED  CF  Ml 

1  T 

J.  •_ 

28.0 

8.0 

35.0 

TO .  0 

CENTER  BELOW 

F.S,  H 

INIMUM= 

719  FOR 

THE 

CIRCLE  OF 

CENTER  <  : 

'7,0,  18.0) 

» 

R 

KR 

BISHOP  MODIFIED  AND/0 

R  ORDINARY 

METHOD  OF 

SLICES 

it- 

1 

perma-f 

-ost  depth 

7  f eat , 

10  feet  from  c 

*~n  1  1  r',E* 

NUMBER 

TANGENT 

RADIUS 

> 

CENTER  ( 

V)  CENTER 

”8 (BISHOP)  FS 

'QMS)  dFS'.BS?  » 

1 

30 . 0 

1 0  0 

35 . 0 

20 . 0 

CENTER  BELOW 

INTERPOLATED  CF.  B 

'“i 

3C' .  0 

1 0 .  j 

7. 1 .  O 

20 . 0 

CENTER  BELOW 

INTERPOLATED  CR,  |<jji 

.616  .000 

V 

30.  0 

1 4 . 0 

35.0 

16.0 

.  668 

4 

30 .  C> 

1 0 . 0 

39.0 

20 . 0 

CF.NTER  BELOW 

INTERPOLATED  CR.  ,v 
INTERPOLATED  CR'  jV 
.  670  .  000 

5 

6 

30 . 0 

30 . 0 

6.0 

14.0 

35 . 0 

33.  0 

24 . 0 
16.0 

CENTER  BELOW 
,779 

30 .  0 

16.0 

35 , 0 

1 4 . 0 

.  682 

, 642  . 000  M 

8 

30 . 0 

14.0 

37.0 

1 6 . 0 

.  646 

.  600  .  000  F? 

9 

30 . 0 

12.0 

35.0 

18.0 

.  667 

.  600  .  000 

10 

30 . 0 

16.0 

37 . 0 

14.0 

.  665 

.  628  . 000  & 

11 

30.  0 

14.0 

39.0 

1 6 . 0 

.668 

,'16  .000  K 

12 

30 .  0 

12.0 

37.0 

1 8 . 0 

.636 

.  580  .  000  kp 

13 

30 . 0 

12.0 

35.0 

18.0 

.  663 

.  600  ,  000  ij 

14 

30.  0 

12.0 

39.0 

18.0 

.  663 

. 600  . 000  " 

15 

30 . 0 

10.0 

37 . 0 

20 . 0 

CENTER  BELOW 

INTERPOLATED  CRi- 

16 

30 .  0 

14.0 

35.  0 

16.0 

.  668 

.616  . 000  w 

17 

30 . 0 

14.0 

39.0 

16.0 

.  668 

.616  .000  ;<■ 

18 

30 . 0 

10.0 

39.0 

20 . 0 

CENTER  BELOW 

INTERPOLATED  CRf  o 

19 

30.0 

1 0 . 0 

35.0 

20 . 0 

CENTER  BELOW 

INTERPOLATED  CR;  j 

F.S.  MINIMUM— 

.  636  FOR 

THE 

CIRCLE  OF 

CENTER  (  77.0,  18.0) 

8 

BISHOP  MODIFIED 

AND/OR  ORDINARY  METHOD  OF 

SLICES 

V 

.*• 

perma-f 

^ast  depth 

7  -feet, 

10  -f 

eet  -from  centerline 

i 

NUMBEF 

TANGENT 

RADIUS 

(/' 

CENTER 

Y )  CENTER 

FS  (BISHOP,'  FS 1  QMS  •  dFS(3S;  £ 

1 

35.0 

15.0 

35 . 0 

20 . 0 

CEM i ER  BELOW 

INTERPOLATED  CRY  .9 

”5 . 0 

15.0 

3  i. .  0 

20 . 

pirMTPP  OJTI  nj.: 

INTERPOLATED  CR;  £ 

“7 

35 . 0 

1 9 . 0 

35-  0 

16-0 

-  60 1 

.  56 6  .  000  | 

4 

35 . 0 

!  5 . 0- 

39,0 

n,*; 

C E N  TE F;  ?  E L  0 ’/ ; 

INTERPOLATED  CF  ~ 

sr 

35.0 

1 1 . 0 

35 . 0 

3  4  -  0 

CENTER  BELOW 

I  ^ITIlFFQL  £TED  Cf  *• 

/. 

*  j 

•33.  0 

19.0 

■,  i 

1  6 . 0 

.  628 

^  '--'OCV  t  ,v 

35 .  0 

2 : .  o 

35.0 

1  '5 . 0 

7 

.37  2  . 002  |* 

o 

u 

35 . 0 

1 9 . 

77, 

1 

J.  '~J  A  J 

.  597 

559  ,  00-*'  9 

9 

35.  0 

1T.0 

35.0 

18.0 

-  oOTj 

.De7  .000 

10 

“rrrr 

-  -J  - 

2 1 . 0 

37 . 0 

1 4 . 0 

,.  596 

.365  .000  I 

1 1 

35  -  0 

19.0) 

79.0 

i  6 . 0 

.  oO  I 

J 

--rr 

1 

37 ,  C 

I  o ,  .‘i 

,  795 

.  55' v  .  5 

J.  -■ 

35 . 0 

2 1 . 0 

35.0 

i  4 . 0 

603 

.372  .000  S 

i  4 

35 . 

2  L  ■  0 

79.0 

1  t  A  0 

,  *0  3 

-5' 72  •  000  V 

15 

35  -  0 

1 7 . 0 

3C ,  0 

t  <z  r, 

,  •  >05 

.36"  .000  | 

16 

•73 .  O' 

1 7 . 0 

75 .  *’• 

18.0 

.  605 

. 767  . OOP  I 

~ .  8 .  i 

UN  TMUM“ 

.5*3  FOF 

T  k” 

CIRCLE  OF 

CENTER  ( 

7 7 .  .  l.i.O- 

*> 

> 

(O  PS-cT  &V=UDV>0  ^oTTona  of  EKASAHtNAS  > 

L. 


Cfc\T\CfcV  fWXOft.  OP  S>M="£TX 


’w.rjirui 


Program  5TABR  —  Version  2,34  IMS-DOS) 


BISHOP  MODIFIED  AND /OP  ORDINARY  METHOD  OF  SLICES 
oerma-f rost  death  1  -feet,  JO  feet  fro m  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  0 
NUMBER  OF  DEPTH  LIMITING  TANGENTS  3 
NUMBER  OF  VERTICAL  SECTIONS  12 
NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  PORE  PRESSURE  LINES  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF  2.0 

ALL  CIRCLES  TANGENT  TO  DEPTH,  28. 0,  30.0,  35.0, 

GEUMhi RY 


SECTIONS 

.  0 

1 0 . 0 

13.  0 

15.0 

20.0 

32.0 

42.0 

64.0 

69.0 

71.0 

74. 

0 

T.  CF:A.CKS 

18.0 

18.0 

18.0 

18.0 

18.0 

18.0 

25 . 0 

25  • 

25 . 0 

3C  - 

—  w  *  V 

75 , 

(■< 

W  IN  CRACK 

18.0 

18.0 

18.  0 

18.0 

18.0 

13.0 

25.0 

75  -  0 

25.  0 

25 . 0 

25. 

0 

BOUNDARY  1 

1S.0 

18.0 

18.0 

18.0 

18,0 

18.0 

25.0 

25.0 

25 . 0 

25.0 

25 „ 

BOUNDARY  2 

25.  0 

25 . 0 

25 . 0 

25.0 

25.  0 

25.  0 

25.0 

25.0 

25,0 

25.0 

nc 

ij 

BOUNDARY  3 

35.0 

35 . 0 

25.0 

25 . 0 

25 . 0 

25.0 

'Eicj  t~f 

25.0 

25.0 

OCT  A 
•- J  * 

’.'W  t 

0 

BOUNDARY  4 

35 . 0 

35 . 0 

28.0 

30 . 0 

35 . 0 

35.  0 

7  s!  i"i 

35.0 

30 . 0 

28.0 

35. 

0 

BOUNDARY  5 

35.0 

35 . 0 

35 . 0 

35.0 

35 . 0 

35.0 

75 . 0 

35 . 0 

35.0 

%!•  -  *J 

— ’cr 

v 

REINFORCING  FORCE  DATA  AT  1  LEVEL <s> 

Y-  25.00  NO.  OF  FORGO  POL-ii  3- 


i  •  •— n 

FORCE 


.  o 

5000 . 0 
5000 . 0 


SOIL  f: POVERTIES 


layer 

i 


COHESION  FRICTION  ANGLE 
'I*  *y,r  t  '  < 

.  '.*00 . 0  40 .  o 

i  .  0  . 

i  000.0  40.0 


INSIT” 


B T SHOP  HOD  'FT ED  AND  'OF  OF  D  fNAP  <  METHOD  >F 
per  mat  r  "r  l  <  i  p  p  L ! ;  *7  f  eec  ,  i  <">  feet  *  rem  :  renter  1  e 
NUMBER  TAMOEY-  RADIUS  CENTER  CENTER 


• ::  i  :hop 


i 

70.0 

12.0 

35 . 0 

18.0 

1  -  305 

1 . 243  . 643  W 

2 

7''.' .  0 

1 2  *  0 

3 1 . 0 

iu  .  0 

1 . 842 

1 . 699NEI5 .  RES  1ST.’"  | 

-r 

30 . 0 

16.  0 

35 . 0 

14.0 

1.317 

i . 277  . 635  ^ 

4 

.70.0 

12.0 

39 . 0 

18.0 

1 . 305 

i  •  24  »•  *  'Id  4  N 

cr 

30.  0 

8.0 

35.0 

22.0 

CENTER 

BELOW  INTERPOLATED  CR£  ? 

6 

30.  0 

12.0 

33.0 

18.0 

1 .  -'74 

1 , 337  .717  j; 

7 

30 . 0 

10.0 

35 . 0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CRf  j- 

8 

30.  0 

14.0 

33 . 0 

16.0 

1  £ 

4  •  • 

1.364  .694  £ 

Q 

30.  0 

1 4 . 0 

3  •? .  A 

16.0 

I. .  268 

1.222  .621  1 

10 

30 . 0 

1  •  j .  0 

37.0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CR.  \ 

l : 

70 . 0 

10.0 

33 . 0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CR:'  ?■ 

12 

30' .  0 

16.  0 

37 . 0 

14.0 

1 . 286 

1,249  .621  ? 

13 

30 . 0 

14.0 

39 . 0 

16.0 

i .  30a 

1 . 254  . 633  [ 

14 

30.  0 

1 6 . 0 

33 .  'I1 

14.0 

1.317 

1.277  . 635  £ 

15 

30 .  O' 

16.0 

39 . 0 

1 4 .  ■:> 

1.317 

1.277  . 635  1 

1 6 

30 . 0 

12.0 

79 . 0 

18.0 

1 . 305 

1 . 243  . 643  t 

•  / 

30 . 0 

12.0 

35.0 

1  O  A 
x  •-/  • 

i  .  305 

S .  Mr 

M  T  M,JM—  i 

^/.p  inp 

THE  CIRCLE  OF 

CENTER  ' 

_■  '  a  ’-!  ■  1 

U  ) 

t 

3 1  SHOP  M 

3D  IF  I ED 

AMD /OP  ORDINARY 

METHOD  OF 

SL ICES 

i 

t 

a 

r  ro-a-r  r 

c--.it  clench 

2  t  r  7 

1 0  r9et  frs:n  •: 

2r ter ; ; no 

* 

MBEF 

TANGENT 

RADIUS 

1 X  >  CENTER 

Y;  CENTER 

- -3  ••  B I 3-:0R  •  ~E  •  QMS;  cF3(3S-  \ 

i 

T3 .  0 

>  7 . 0 

,  %  i  i 

18.0 

.  83“ 

. 332  . 22c  f 

-•■cr  t*( 

1 " .  0 

i  o .  ■:> 

■  r 

2 

^  1  i*. 

~“TT 

i  .T. 

1  ^ 

.  -  i  9  .  T  1  * 

4 

it 

« 7  r. 

7^ .  v 

:  o  *“• 

38“ 

Tf*-  -per  t 

-*r  t\ 

*  *■*’ 

— —  - 

'I'- 

rr-r- 

•  - 

*-  -  *  • 

‘5 

.  — ^  . 

i  ”  % 

3 

r  #*rr 

*3 - 5  .  P 

7 

T3 .  0 

•  **•*  |  « 

'  T. 

i  ::  -  * ! 

.  C  * 

a -a  1 

‘  '  1  -  t 

15 

T5.  0 

*  “7  1*1 

i  >  <  > 

4  1  • 

T  —  r* 

.  3°7  - 2vc 

<  «"« 

1.  j  « 

T  - .  \  • 

r*  * " 

<  ■  “  r* 

~J. F>_ ';W  F  'iTF“'r:0i.Ai  i'ED  CR 

i  v 

—  cr 

1  "*  .  V 

.  '*;  <  1 

l  n  ,\ 

.252  .  ST  • 

i  : 

:  Q  *‘t 

.  ; 

!  ^ 

,  *■  i 

15.0 

“'l*l  *1 

CE.’iT'TF 

BELOW  1  i  ITERPOLATED  C~ 

L  7 

75 . 

1 9 . 0 

*’3- 

1  t  t^l 

. -  70 ~ 

Program  STftBR  —  Version  2.84  (MS-DOS) 

BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  -feet,  10  test  -from  centerline 
CONTROL  DATA 

NL'MBER  OF  SPECIFIED  CENTERS  0 

NUMBER  OF  DEPTH  LIMITING  TANGENTS  1 

NUMBER  OF  VERTICAL  SECTIONS  12 

NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 

NUMBER  OF  PORE  PRESSURE  LINES  0 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 

SEISMIC  COEFFICIENT  SI, 52  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 

SEARCH  STARTS  AT  CENTER  <  37.0,  22.0)  WITH  -INAL  GRID  OF 


GEOMETRY 

-w* 

SECTIONS 

•*' 

10.  0 

13.0 

15.0 

T.  CRACKS 

18.0 

18.0 

18.0 

18.0 

W  IN  CRACK 

18.0 

18.0 

18.0 

18.0 

BOUNDARY  1 

18.0 

F** 

CO 

o 

18.0 

18.  0 

BOUNDARY  2 

25.0 

25.0 

25.0 

25 . 0 

BOUNDARY  3 

35.0 

35.0 

25 . 0 

25.  '! 

BOUNDARY  4 

35.0 

35 . 0 

28.0 

30. 0 

BOUNDARY  5 

35 . 0 

35.0 

35.0 

35 . 0 

REINFORCING  FORCE  DATA  AT 

1  LEVEL 

; 

23 , 0 
25  .  o 


25.0 
25 . 0 
25.0 
25 . 0 
25.0 
28.0 
35 . 0 


74.0  * 


r~,  t—  Jl 

2d  .  ! 

25.0  jf 
25.0  <; 

o5 .  v  i 
35 .  0 


25 . 00 
X 

42-  0 
79.0 
r;  i*' 


NO.  OF  FORCE  --O' 
FORCE 
.  0 

7000, 0 
7 000 . 0 


SOIL  PR OPERATES 

LAVER 

1 


COHESION  FRICTION  ANGLE 


-  -J  -  V 

u.{  ■  «  i 


ICTTV 


i  000 . 0 
i  0*"‘  •  0 
1000. 0 


B  r  SHOP  MODIFIED  AND 'CP  OR  ~  i  MARY  f-ET'-'OD  OF  21  ICES 
oermafrost  depth  7  reef  ,  1  *1*  feet  ft  om  “enter- 1  i  v- 
NUMBER  TANGENT  RADIUS  •'  < ;  CENTER  (Y /  CENTER  FS  B  7  "HOP 


C-J  K> 


Program  STA8F,  —  Version  2.84  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  -Feet,  10  feet  -from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  0 
NUMBER  OF  DEPTH  LIMITING  TANGENTS  1 
NUMBER  OF  VERTICAL  SECTIONS  12 
NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  PORE  PRESSURE  LINES  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  3i,S2  =  .00,  .00 

UNI'i  WEIGHT  OF  WATER  -  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF 

ALL  CIRCLES  TANGENT  TO  DEPTH,  35.0, 

GEOMETRY 


SECTIONS 

.  0 

10.  0 

1 3 . 0 

15.0 

20.  0 

32 . 0 

42.0 

64 . 0 

69.  0 

7 1 . 0 

74.0 

T.  CRACKS 

18.0 

18.0 

18.0 

L  8 . 0 

18,0 

18.0 

25.0 

25.0 

25 . 0 

25 . 0 

25,  0 

W  IN  CRACf: 

18.0 

18.0 

18.0 

18.0 

18.0 

13.0 

25 . 0 

25 . 0 

25. 0 

25.0 

25,0 

BOUNDARY  1 

18.0 

J  8.  0 

18.0 

18.0 

18.0 

18.0 

25.  0 

25 . 0 

23 . 0 

25.0 

25 .  O 

BOUNDARY  2 

25 . 0 

25 . 0 

25.  0 

25.0 

25.  0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

BOUNDARY  3 

35.0 

35 . 0 

25.0 

25 . 0 

25.0 

25 . 0 

25  „0 

25 . 0 

25 . 0 

23.0 

35,  0 

BOUNDAF;Y  4 

35 . 0 

35.  0 

28 . 0 

30,0 

35 . 0 

35.0 

35.  0 

35 . 0 

30. 0 

28. 0 

35 . 0 

BOUNDARY  5 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35.0 

35 . 0 

35 . 0 

35.  0 

3j,  <"■ 

35 . 0 

35 . 0 

REINFORCING  FORCE  DATA  AT  1  LEVEL (3) 


25.00 

X 

42.0 
39 . 0 


NO.  OF  FORCE  POINTS* 
FORCE 


,  0 

9000. O 
9Ou0 -  A 


BOIL  PROPERTIES 


avl-:r 

1 

CQHEF 

, 

I  ON 

FRICTION  ANGLE 

33 . 0 

DENSITY 
J  35.0 

z 

1 OOO . 

'3-0 . 0 

1  20 . 0 

i  o(  > .  < 

.  0 

i  or. ,  0 

A 

1000.  ( 

1 

40 . 0 

1 20 . 0 

B'BHDP  MODIFIED  AimD/QR  ORDINARY  METHOD  OF  SLICES 
perm*-*  !"C£,L  depth  7  feet.,  1 0  feet  from  cent'-r  1  inv 
NUMBER  TANGENT  F.iDRJS  'X-  CENTER  <\)  CENTER  PS- PI SHOP; 


FS t QMS) 


dFS'DS' 


TOnGViLzx  smn’jfn&fjz  / : 

i 

.  ■ .' 

1  ■  >  Q 

37 . 

35 . 0 

13.0 

-*■ 

35.  0 

17.0 

-•  /  » 

4 

35.0 

i  3 . 0 

41 

5 

35 .  0 

9 . 0 

“^7 

6 

35 . 0 

1 7 , 0 

“r(T 

7 

35.  0 

19.0 

3  / 

S 

35 .  0 

5  7.0 

39 

9 

35.  0 

15.0 

-y? 

10 

35 , 0 

1 9 . 0 

7Kr. 

1  t 

35 . 0 

19.0 

79 

J.5-?  J^^TjaWTnR-TLTOJ  -.^’l 


0 


program  STABR  --  Version  2. 84  (MS-DOS) 


BISHOP  MODIFIED  AMD/UR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  -feet,  JO  -feet  from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


O 


t  — 

5 

0 


.  00 , 


,  00 


SEISMIC  COEFFICIENT  S1,S2 
UNIT  WEIGHT  OF  WATER  -  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  <  37.0,  22.0)  WITH  FINAL  GRID  OF  2.0 

ALL  CIRCLES  TANGENT  TO  DEPTH.  35.0, 


GEOMETRY 


SECTIONS 

.  0 

10.0 

13.  0 

15.0 

20 . 0 

32.0 

42 . 0 

64 . 0 

69 . 0 

7 1 . 0 

74.  0 

T.  CRACKS 

18.0 

18.0 

18.0 

18.0 

18.0 

18,0 

25.  0 

25.0 

>,io .  0 

25 . 0 

25.  0 

W  IN  CRACK 

18.0 

18.  0 

18.0 

18.  0 

18.0 

18.0 

25,  0 

25.0 

25 . 0 

25.0 

BOUNDARY  1 

13.0 

18.0 

18.0 

18.0 

18.0 

18.0 

25.0 

25. 0 

25.0 

25 , 0 

r. 

■l_  W  0  - 

BOUNDARY  2 

25.  0 

25  .  0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

•1 GT 

25 . 0 

25 . 0 

BOUNDARY  3 

35. 0 

35 . 0 

25 . 0 

25 . 0 

25 . 0 

25 . 0 

25.0 

25.0 

25.  0 

25.0 

35.  0 

BOUNDARY  4 

35.0 

35 . 0 

28.0 

30 . 0 

35 . 0 

35.  0 

35.0 

35 . 0 

30.  0 

28 . 0 

35.0 

BOUNDARY  5 

35.  0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35 , 0 

35.  0 

35. 0 

35.0 

33 . 0 

35 .  -I' 

REINFORCING  FORCE  DATA  AT  1  LEVEL (s) 


Y* 


23 ..  00 


X 

42. 0 
39 . 0 


NO.  OF  FORCE  POINTS 
FORCE 
.  0 

ROOO . O 
8000 . 0 


=  4 


0 


SOIL  PROPERTIES 


LAYER 

1 


4 


COHESION  FRICTION  ANGLE 
33 . 0 

1 000 . 0  40 . 0 

100.0  .0 

i. . 


l  i 


f 

'> 


L  21"' .  0 


looo. o 

BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
perma-frost  depth  7  -feet,  J  o  feet  from  r  enter  1  j  ::e 
NUMBER  rAMGEMT  RADIUS  -  X )  CENTER  (Y;  CENTER  F“<BTSH:jt: 


S ■ QMS ) 


ClFS  •  B,: 


i 

■Zu .  <J 

1 3  m  (_i 

37  *  i  > 

22.0 

CENTER 

BELOW  INTERPOLATED  CF  KJ 

35,  0 

1 3 . 0 

33.0 

22.0 

CENTER 

BELOW  INTERPOLATED  CP  jCM 

35.0 

17.0 

37 . 0 

18.0 

1 . 043 

1 . 007  , 44*  B 

tX 

35  n  0 

1  3 . 0 

9 1 ,  '"J 

22..  0 

CENTER 

BELOW  INTERPOLATED  CR  SS 

rrt 

35,  0 

9.0 

37 . 0 

28 . 0 

CENTER 

BELOW  RIGHT  INTERSECT  f* 

6 

35 . 0 

17,0 

35..  0 

18.0 

1 « 060 

J. .  023  .  45c  $ 

7 

35.  0 

1 9 . 0 

37 . 0 

1 6 . 0 

1 . 072 

1 . 038  , 47*  m 

8 

35 . 0 

J  7 . 0 

39,0 

.'8.0 

i  ,  060 

j . 023  „ 456  W 

9 

35 . 0 

1 5 . 0 

37 . 0 

20.  0 

CENTER 

BELOW  INTERPOLATED  2F  gj 

10 

35 . 0 

19,0 

35 . 0 

1  6 . 0 

1 . 086 

1,051  ,485 

1 1 

35.0 

19.0 

39. 0 

1  £  „  0 

1 , 086 

1,051  .  -85  ($ 

•*  '**) 

L  4.. 

35.  0 

15,0 

39. 0 

20 .  0 

CENTER 

BELOW  INTERPOLATED  CR 

i  3 

35.  0 

15.0 

35 . 0 

zo  *  0 

CENTER 

BELOW  INTERPOLATED  IR.  || 

s. 

M I N I MUM- 

1.043  FOR  THE 

CIRCLE  OF 

CENTER  ( 

37 . 0 ,,  j  8 

.  o )  m 

^  ", 

rf. 

£ 


u>: 


V  . 
<•*, 
v, 

y* 

i 


^v.v.y.y. 


’/csrss^s*  .j,j~*i7^s^x^;;^or«A 


F'roqram  STABR  —  Version  2.84  (MS-DOS)' 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  Oepth  7  feet.  10  feet  from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


SEISMIC  COEFFICIENT  Si,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF 

ALL  CIRCLES  TANGENT  TO  DEPTH,  28. 0,  30.0,  35.0, 

GEOMETRY 


SECTIONS 

.  0 

1 0 . 0 

1 3 . 0 

15.0 

20. 0 

32.  :• 

42.0 

64 . 0 

69.  0 

71.0 

74.0 

T.  CRACKS 

18.0 

18.0 

18.  C 

18.0 

IS.  0 

18.0 

25.0 

25.0 

25 . 0 

25.0 

25.  0 

W  JN  CRACK 

18.0 

18.  0 

1 8 . 0 

18.0 

18.0 

18.0 

25 . 0 

25 . 0 

25.0 

25 . 0 

25..  0 

BOUNDARY  1 

IS.  0 

3  3.0 

18.0 

18.0 

18.0 

18.0 

25 . 0 

25 . 0 

25. 0 

25.0 

25 . 0 

BOUNDARY  2 

25.  0 

25.  0 

25.0 

25.0 

25 . 0 

25 . 0 

25.0 

25. 0 

25.  0 

25 . 0 

~5.  U 

BOUNDARY  3 

35 . 0 

35.  0 

25.0 

25.  0 

25.0 

25 . 0 

25 . 0 

25 . 0 

25.0 

25 . 0 

35.  0 

BOUNDARY  4 

35.0 

35.  0 

28.  0 

30.0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

30 . 0 

28.0 

35 . 0 

BOUNDARY  5 

35.  0 

35.0 

35 . 0 

35.0 

35. 0 

35 , 0 

35 . 0 

35.0 

35 . 0 

35.0 

35.0 

REINFORCING  FORCE  DATA  AT  2  LEVEL (s) 


C„-  .  ‘-’V 

X 

37.7 
34 . 0 
3 « 0 
-  0 


NO.  OF  FORCE  POINTS- 
FORCE 

.  0 

3000 . 0 
5000. 0 

.  0 


NO.  OF  FORCE  POINTS* 
FORCE 


3000 .  0 
5000. 0 

.  0 


on  r  i  osrnirr.T  r 


LAN  ER 


COHESION 

.  0 

i  000 .  0 

1 0'j , 1 


FRICTION  AMFLr 
r11*  r  i 


Dt’-iSi  TV 


tx 

1  ooo . 

0 

AO.O 

1 20 . 0 

fc 

i 

fc 

B I SHOr-'  MODIFIED  f 

iND/OR  ORDINARY  METHOD  OF 

SLICES 

permaf  .* 

-  ost  depth 

7  -feet, 

10  f 

eefc  from  centerline 

\ 

i 

i 

NUMBER 

TANGEN1 

RADIUS 

v  X ) 

CENTER 

( Y  i  CENTER 

i-S  (BISHOP.'  FS  '.QMS;  LFS'.BS 

1 

28.0 

6 . 0 

37.0 

22 « 0 

CENTER 

BELOW  INTERPOLATED  CF  jj 

T1 

28 . 0 

6  >  0 

33 . 0 

22 . 0 

CENTER 

BELOW  INTERPOLATED  CP  | 

“T‘ 

20.0 

i  0 .  'I1 

37.0 

i  3 . 0 

•  -!  — \  *_  c-; 

2.  IG^  L.34-.  | 

& 

28 . 0 

6  0 

4  i  ,  0 

22.  0 

CENTER 

BELOW  INrERPOLATED  CP  i 

5 

28 . 0 

2  -  0 

37 . 0 

26 . 0 

CIRCLE 

OUTSIDE  SLOPE  jj 

6 

28.  0 

1 0 . 0 

35 . 0 

18.0 

2.404 

2. 3]. 2  1 . 632  S 

7 

28.0 

12.0 

37 . 0 

1 6 . 0 

2.  307 

2.240  1.567  j 

8 

28.0 

1  0 . 0 

39 . 0 

1  8 . 0 

2  n  *1 0  4 

2.312  1.632  ! 

9 

28.0 

8.0 

37 . 0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CR  l 

10 

28.0 

12.0 

35 . 0 

1 6 . 0 

2.41 0 

2.341  1.631  j 

1 1 

28.0 

L  2 . 0 

39 . 0 

16.0 

2 .  41 0 

2 .341  1 . 63 i  j 

12 

28.0 

8.0 

39 . 0 

20.0 

CENTER 

BELOW  INTERPOLATED  CP:  [ 

13 

28.  0 

8.0 

35 . 0 

20.0 

CENTER 

BELOW  INTERPOLATED  CRY ! 

F .  S .  M 

INI MUM  =  2 

.265  FOR 

THE 

CIRCLE 

OF  CENTER  < 

37 . 0 ,  18. 

0  >  r 

BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF 

SLICES 

permaf 

rash,  depth 

7  feet, 

10  f 

eet  from  centerline 

__  _ _ _ 

NUMB '^n' 

i ANGENT 

RADIUS 

<  X ) 

CEN  i  tip 

( Y )  CENTER 

FS  ( E I SHOP '  FS ( QMS ;  dFS ( BS ;  - 

<k 

1 

30 . 0 

10.0 

35.0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CRT ! 

*~t 

30 . 0 

1 0 . 0 

31.0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CRT- j 

3 

70.0 

1 4 . 0 

35.  0 

1 6 » 0 

1 , 73 1 

j. .  680  1 , 063 

4 

30 . 0 

1 0 . 0 

39 . 0 

20.0 

CENTER 

BELOW  INTERPOLATED  CP:-  1 

5 

30. 0 

6.0 

35 . 0 

24.  0 

CENTER 

BELOW  INTERPOLATED  CRT 

6 

30.  0 

14.0 

33 . 0 

16.0 

1 . 895 

1.82o  1.156  1 

7 

30 . 0 

1 6 . 0 

35  n  0 

1 4 .  <'• 

1 . 779 

1 . 779  i . 09  7 

0 

30 . 0 

1 4 . 0 

37  -  0 

16.0 

1 . 682 

1.636  1.035 

9 

30.  0 

12.0 

35 . 0 

1 8 . 0 

1 . 672 

1.610  1.010 

10 

30 . 0 

12.0 

33.0 

1 8 . 0 

1 . 884 

L . 79”  1.127 

1. 1 

30 . 0 

12.0 

37 . 0 

18.0 

1 . 612 

1 . 556  . 976 

1  o 

30 . 0 

1 0 . 0 

35 . 0 

20.0 

CENTER 

BELOW  INTERPOLATED  CPf 

1  _T 

T'-‘ .  0 

14..  ft 

•77.0 

1 6 . 0 

1 . 682 

1  >  67*6  '  »  0“2*“ 

14 

3  r,  m  r, 

12.0 

39.0 

18.0 

1 . 672 

1  „  6  i 1  f  1  „  0  i  0 

15 

30 . 0 

1 0 . 0 

37 . 0 

20.0 

CENTER 

BELOW  INTERPOLATED  CR: 

16 

30 . 0 

1 4 . 0 

35 . 0 

1  o .  0 

i .  ”31 

1 . 6S0  i . 06 : 

17 

30 , 0 

J4,r. 

39.0 

16.0 

1.731 

1 . 65'"1  1. 063 

18 

30 .  l' 

1 0 . 0 

3  .  0 

20 . 0 

CENTER 

BELOW  INTERPOLATED  CP 

19 

Z* j  r' 

:  0  „  0 

35 . 

20 , 0 

CENTER 

BELOW  INTERPOLATED  CP 

F.S.,  !- 

INI MUM-  1 

.112  FOR 

IKE 

CIRCLE 

OF  CENTER  < 

~'r  ~7  m  -I 

-  ■  p  “  7  - 

,  0 ) 

u^,Lriu/io-T«r.  ’■*.*/■»  imrjjv  *jvajwxwxj*n«n^jSWT  t\tl 


Program  3TADR  —  Version  2.8*1-  .'MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
jsrM-fr ost  dep th  7  Teet,  10  -feet  -frcsm  :enteri  me 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  O 
NUMBER  OF  DEPTH  LIMITING  TANGENTS  3 
NUMBER  OF  VERTICAL  SECTION'S  1 2 
NUMBER  QF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  POPE  PRESSURE  LINF.S  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 


SEARCH  STARTS  AT 

CENTER  (  3 

7.0, 

22.  0) 

WITH 

FINAL 

GRID  QF 

2.0 

ALL  CIRCLES  TANGENT  TO 

DEPTH 

,  28. 

0,  3 

0.  0, 

35.  0, 

GEOMETRY 

SECTIONS 

.  0 

10.  0 

13.0 

15.0 

20 . 0 

32.0 

42.0 

64 . 0 

69 . 0 

71.0 

T.  CRACKS 

18.0 

18.  0 

18.0 

18,0 

18.0 

13.0 

25.0 

25.0 

25,0 

25 . 0 

W  IN  CRACK 

13.0 

18.0 

18.  0 

18.0 

1 8 . 0 

18.0 

25.0 

25.  0 

25.  0 

25.0 

BOUNDARY  1 

IS.  0 

18.0 

18.0 

13.0 

18.0 

18,0 

25,  0 

25.0 

25.0 

25.0 

BOUNDARY  2 

25.0 

25.  0 

25.  0 

25 . 0 

25.0 

25,  0 

25.0 

25.  0 

25 . 0 

25 . 0 

BOUNDARY  3 

35.0 

35.  0 

25 . 0 

25.0 

25.0 

25 . 0 

25 . 0 

25 . 0 

25.0 

25 , 0 

BOUNDARY  4 

35.  0 

35.  0 

28.0 

30.  0 

35 . 0 

35.0 

35. 0 

35 . 0 

30.  0 

28.0 

BOUNDARY  5 

35.0 

35 . 0 

35.0 

35.0 

35 . 0 

35.  0 

35.0 

35.0 

35. 0 

_J  M  O 

REINFORCING  FORCE  DATA  AT  3  LEVEL'S) 


20.  00 
X 

:.4„8 
:•  1 . 0 


MO.  OF  FORCE  POINTS” 
FORCE 
.  0 

1666. 0 
1 666 . 

.  0 


22.  00 

X 

37, 7 
34 . 0 


NO.  OF  FORCE  POINTS1 
FORCE 


5O00. 
50»''0 . 0 

.  0 


23.00 

X 

a  Zl  •  0 
39,0 
3.  O 

("i 


NO.  OF  FORCE  POINTS- 
FORCE 


5000. O 
5000 . 0 

.  0 


c/i  Li  i.n  La  !_a 


tes.  ■^jrJ^35{JESjHajix3^^oc.<jE?«^Ka^,T^t^wz3c^«,r^2CT^<^xjTCT^rss?cy.7rsvJ 


SOIL  PROPERTIES 


LAYER 

COHESION  FRICTION  ANGLE 

DENSITY 

1 

o 

35.0 

1 35 . 0 

o 

1 000 

0 

40 . 0 

1 20 .  ■:> 

T 

1  C'O . 

0 

.  0 

3  05 , 0 

4 

1000. 

o 

40.0 

i  20 . 0 

BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF 

SLICES 

permafrost  depth 

7  feet, 

10  -feet  -from  centerline 

NUMBER 

TANGENT 

RADIUS 

( X )  CENTER 

fY)  CENTER 

FS  ( B  J  SHOP  -  FS  <■  0M5 )  dFS  <  BS  > 

i 

2S.0 

6. 0 

37.0 

22 . 0 

CENTER  BELOW  INTERPOLATED  CF- 

n 

28.0 

6.  0 

33.0 

22.0 

CENTER  BELOW  INTERPOLATED  Cci 

"31 

28.0 

1 0 . 0 

37 . 0 

18.0 

2,358  2.283  1.639 

4 

28.0 

6.  0 

41  .0 

22. 0 

CENTER  BELOW  INTERPOLATED  CEr 

5 

28.  0 

2.0 

37.  0 

26 . 0 

CIRCLE  OUTSIDE  SLOPE 

6 

28.0 

i  0 . 0 

35.  0 

18.0 

2.503  2.411  1.771 

7 

28.0 

12.0 

37.0 

16.0 

2.446  2.388  1.706 

S 

28.0 

10.0 

39 . 0 

3  3.0 

2.503  2,411  1,731 

9 

28.0 

3.  0 

37 . 0 

20 . 0 

CENTER  BELOW  INTERPOLATED  CFT 

3  0 

28.0 

12.0 

35 . 0 

3  6 . 0 

2,555  2.486  3,776 

11 

2B .  0 

12.0 

39.0 

3  6 . 0 

2.555  2.486  1.776 

12 

28.0 

8. 0 

39 . 0 

20.0 

CENTER  BELOW  INTERPOLATED  CRY 

13 

28.0 

8.0 

35 . 0 

20.0 

CENTER  BELOW  INTERPOLATED  CF/- 

F .  S .  M 

INTMUM"  2 

.358  FOR 

Tf  *  CIRCLE 

OF  CENTER  < 

37.0,  18.0 i 

BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF 

SLICES 

perma-f 

rast  depth 

7  feet, 

10  feet  from  centerline 

NUMBER 

TANGENT 

RADIUS 

(X)  CENTER 

<Y ■  CENTER 

FSt BISHOP)  FSvOMS)  dFS'BS- 

1 

30 . 0 

1 0 . 0 

35 , 0 

20.0 

CENTER  BELOW  INTERPOLATED  ORf 

30.  0 

1  0 . 0 

31 . 0 

20 . 0 

CENTER  BELOW  INTERPOLATED  CPf 

"7 

30.  0 

14.0 

35.0 

16. 0 

1.826  3 . 774  1.138 

4 

30.0 

1 0 . 0 

39 . 0 

20.0 

CENTER  BELOW  INTERPOLATED  CPf 

0 

30 . 0 

6 . 0 

CT 

24 . 0 

CENTER  BELOW  INTERPOLATED  CRf 

6 

30.  0 

1 4 , 0 

33.  0 

16.0 

1 . 998  1 . 929  L . 259 

~7 

“V.  % 

16.  0 

35 « 0 

1 4 . 0 

1.895  1.854  1.212 

3 

30 . 0 

14.0 

37 . 0 

1 6 . 0 

1.774  1.728  1.127 

9 

30 . 0 

12.0 

35,.  0 

1  8 .  A 

i  ,  774  3. ,  67 1  1 . 07  i 

1  i  l 

30-  0 

12.0 

33 . 0 

13.0 

L.952  ’.365  1,195 

1  1 

3 '  • ..  0 

12,0 

33 

3.8.0 

3  .  67  !  1 ,  {■■  15  1  /'72 

i 

30-  0 

J.  v « •  * 

35 .  f» 

20 . 0 

CENTER  D'EL.OW  i'MrEF.f-'OLATED  2F. 

L  '■'-' 

30 . 0 

J.  4 .  o 

37  -  0 

1 6 . 0 

1 . 774  ,  728  L  .  3.27 

14 

30  -  0 

1 2 . 0 

39.0 

18.  v 

i . ”34  -.271  1.071 

i  b 

30  N  0 

1  o .  o 

Z7 . 0 

20.0 

CEH’EF'  BELOW  INTERPOLATED  CR 

16 

30.  0 

L  >' ,  0 

35.0 

'  6 . 

L . 826  1 . 77-  l . i 58 

1 7 

30  „  0 

»  ."i  ■  » 

39 . 0 

L  Cl*  -  *  ’ 

13 

30-  0 

i  o .  ■:> 

~9  „ 

20 . 0 

CEHTEE  BELOW  INTERPOLATED  CF 

1? 

3»I’  0 

3  0 , 0 

*»  m  1 1 

21 > . 0 

CENTER  BELOW  INTERPOLATED  CP 

F.S.  MINIMUM*  1 

.671  FOR 

THE  CIRCLE 

OF  CEN  uIR  ' 

N 


R I SHOP  MOD I TIED  AMD/OP  ORDINARY  METHOD  OF  SLICES 
psnsafrojt  o^pth  ;  ee'r. ,  10  feet  troin  cent'-’r  1  j  ncv 


-£  *  "'-A  ^ 


Z!«n'.vw«'  rxv*s*-V7*i:~*nrr,ir"-ir*  w.-~*  {.n^—  .ctni^s-csac-  j 


Program  STABR  —  Varsion  2.84  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  feet,  10  feet  from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


o 

12 

5 

0 

0 


SEISMIC  COEFFICIENT  S1,S2 
UNIT  WEIGHT  OF  WATER 


. 00 ,  . 00 
62.40 


SEARCH  IS  BASED  DM  BISHOP  MODIFIED  METHOD 


SEARCH  STARTS  AT  CENTER  (  07.0,  22.0)  WITH  FINAL  GRID  OF  2.0 


ALL  CIRCLES  PASS  THROUGH  THE  POINT  (  42.0. 


i  "i  \ 


GEOMETRY 


SECTIONS 

.  0 

10.  0 

13.0 

o 

in 

20 . 0 

32 . 0 

42.0 

64.0 

6? .  0 

7 1 . 0 

T.  CRAC1 S 

CO 

o 

18.0 

18.0 

18.0 

18.0 

18.0 

25.0 

25.0 

25 . 0 

23 . 0 

W  IN  CRACK 

18.0 

18.0 

18.  0 

18.0 

18.0 

18.0 

25.0 

25.0 

25.  0 

25.0 

BOUNDARY  1 

18.0 

18.0 

18.0 

1S.0 

18.0 

18.0 

25.0 

25.0 

25 . 0 

25 . 0 

BOUNDARY  2 

25.  0 

25.  0 

25.  0 

25.0 

25.0 

25.  0 

25 . 0 

25.  0 

25.0 

25.  0 

BOUNDARY  3 

35.0 

35 . 0 

25.0 

25.0 

25.0 

25.0 

25.0 

25 . 0 

23.0 

7:5 . 0 

BOUNDARY  4 

35.  0 

35 .  0 

28.0 

30 . 0 

35.  0 

35.0 

35.0 

35.  0 

30.0 

28.0 

BOUNDARY  5 

35.0 

35.0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

74.0 


REINFORCING  FORCE  DATA  AT  0  LEVEL (s) 


SOIL  PROPERTIES 


LAYER 

1 


COHESION  FRICTION  ANGLE 


.  v 


1 00< 
10( 
1  ooc 


• .  o 

>.  0 
> .  o 


.-5 .  u 
40 . 0 
.  0 
4’!» .  0 


DENSITY 
1  35 . 

1 20 . 0 
105 . 0 
120.0 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 


NUMBER 


4 

rr 

J 


et  depth 

■£  ri  o  . 

10  feet  from 

c  sn  t.  1  ir.0 

T  ANGENT 

PAD I US 

•'  y )  CENTER 

{ v  ■  r* cm 7 err; 

'  T  c_l  1  1 

~B  ■  B 7  SHOP  •  FC  'CMS'  dFE  • 

DC 
—  '-1 

35.? 

13.9 

37 .  v 

,*t 

CENTER  BELOW  INTERPOLATED 

Cr 

-  O 

15.8 

37 , 0 

, 

CENTER  BELTW  lliTEFF CLAWED 

CF* 

7  7 

17.7 

3  '  •  0 

IT." 

■  J  '  *  .'  1  « 

vO'. 

-■5  • 1 ' 

i  ‘  o 

4  j  -  V 

*'  '“t 

CENTFF  BELOW  I  MTEB ::,OLATED 

pr; 

T  6 «  3 

10.  " 

3~ .  0 

26.  0 

CENTER  BELOW  *  ■  GRT  -NT EPS 

ECT 

76.  4 

1 8.  1 

1 8 .  i 

35.  6 

i r- .  c 

37.0 

i  ::  •  3* 

.  !  1 0 

~3  ^  7 

17-3 

W  1  < 

!  r.  .  " 

rwv 

>:  i  ■  j  vi  t  j  ro  m  ro 

CO  1.0  CO  CO  CO  CO  C.0 


n  n 


i~t\LV\i.'.yU\'  J\'J  1  -J  :rer~*_*^— » r>J3CT.v— ,rr-i.~*-vr»  '-"^-V"*  st»  i~»  v*  •-•»  t~»  \r*  \r*  \n  i-w 


Progr am  3TABF.  —  Version  2.84  (MS-DOS) 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  feet,  10  feet  from  center 1 i re 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS 

NUMBER  OF  DEPTH  LIMITING  TANGENTS 

NUMBER  OF  VERTICAL  SECTIONS 

NUMBER  OF  SOIL  LAYER  BOUNDARIES 

NUMBER  OF  PORE  PRESSURE  LINES 

NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE 


o 
0 
i  n 


0 

0 


SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF 

.ALL  CIRCLES  PASS  THROUGH  THE  POINT  (  42.0, 


!>■ 


>5  .  O  l 


GEOMETRY 


rx 


SECTIONS 

• 

f'l 

10.0 

13.0 

15.0 

20 . 0 

32.0 

42.0 

64. 0 

69 . 0 

71 . 0 

74.0  K' 

T.  CRACKS 

18. 

0 

18.0 

18.0 

18.0 

1 8 . 0 

18.0 

25.0 

25  ■  •  j 

23 . 0 

23.0 

.  S 

,  l*  jm 

W  IN  CRACK 

18. 

o 

18.0 

18.0 

18.0 

18.0 

18.0 

25.0 

25.0 

25.0 

25.0 

25  *  0  p? 

BOUNDARY  1 

18. 

0 

18.0 

18.0 

1.8.0 

18.0 

1 8 .  •  j 

25.0 

f. 

23 . 0 

r, 

ww/  i  J  . 

BOUNDARY  2 

25* 

•I) 

25.0 

25.0 

25.0 

25.0 

25.  0 

25,0 

25,0 

23.  0 

25.0 

23.0 

BOUNDARY  7 

35* 

0 

35.0 

25 . 0 

25 . 0 

25 . 0 

25.0 

OCT 

23.0 

ocr 

L  -•  *  J 

oc;  r. 

35.0  £ 

BOUNDARY  4 

—tr 

0 

35.0 

28.0 

30 . 0 

35 . 0 

35  -  0 

35.0 

35 . 0 

30 .  0 

28. 0 

35.0  & 

BOUNDARY  5 

33. 

o 

35. 0 

35 . 0 

35. 0 

35.0 

35.0 

35.0 

35 . 0 

33. 0 

35 . 0 

35.0  | 

REINFORCING  FOR 

CE  DA 

TA  AT 

1  LEVE 

L  (s) 

25 . 00 

NO.  OF 

FORCE 

POINTS 

=  4 

„  » 

!v 

.v 

v 

r\ 

42.0 
:•? .  0 
3 . 0 

.  0 


FORCE 

.  0 

5000 . 0 
3000 . 0 

.  0 


i 


tV 


SOIL  PROPERTIES 
■LAYER 


& 

1 


COHESION  FRICTION  ANGLE  DENS  IT.' 

~3.v  :  37 ,  *'• 

40 .  r.  1 2  .  0 

a.'r. 

BISHOP  MODIFIED  AND /OF  OFC  ;MAP;  ^ET-OD  OF 

t  -Gp  tn  7  f  ^  i  *r.pF  -f, {?n ii sr  1 1  * ^ 


i  •.*  w . v 
1  vO  .  0 
1000*  0 


: 1  Ti , 1  • 
1  20 . 0 


7V 

•V 

v; 

? 


->  i _ 


NUMBER  TANGEfi' 


r  ■  -\  \r-p 


i."JZ  kX^JT: 


-c  .vr*  y;xi?r.vvLi 


1 

35 . 9 

13,9 

37.0 

22.0 

37. 8 

15.8 

33-  0 

22.  0 

■7 

35.7 

17.7 

37.0 

18.0 

4 

35 . 0 

13.0 

41.0 

22 .  *“• 

5 

36. 3 

i  0 , 3 

37 . 0 

26. 

6 

36  •  4 

18.4 

rr  *. 

18.0 

7 

35.  6 

19,6 

37.0 

1 6  „  0 

8 

rc  ■“*■ 

.•%i  •  -* 

17.3 

39 . 0 

5.0 

O 

35.8 

15.8 

37 . 0 

_0.  0 

1 0 

37.2 

19.2 

33 . 0 

18.0 

1 1 

36.2 

20 . 2 

35. 0 

16.0 

i  2 
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18.0 
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16.0 
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33.0 

20 . 0 

16 

38 . 0 

22.0 

7 1  (“• 

16.0 

17 

36 . 2 

20 . 2 

35. 0 

16.0 

13 
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16.6 
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20.0 
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I#  Program  STABR  -■■  Version  2.34  (MS-DOS) 


BISHOP  MODIFIED  AND. OR  ORDINARY  METHOD  OF  SLICES 
permafrost  depth  7  feet,  10  feet  -from  centerline 


CONTROL  DATA 

NUMBER  OF  SPECIFIED  CENTERS  0 
NUMBER  OF  DEFTH  LIMITING  TANGENTS  0 
NUMBER  OF  VERTICAL  SECTIONS  12 
NUMBER  OF  SOIL  LAYER  BOUNDARIES  5 
NUMBER  OF  PORE  PRESSURE  LINES  0 
NUMBER  OF  POINTS  DEFINING  COHESION  PROFILE  0 


SEISMIC  COEFFICIENT  S1,S2  =  .0I>,  .00 

UNIT  WEIGHT  OF  WATER  =  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF 

ALL  CIRCLES  PASS  THROUGH  THE  POINT  (  42.0,  30.0) 

GEOMETRY 


SECTIONS 

.  0 

i  0 . 0 

13.0 

15.0 

20 . 0 

32.0 

42 . 0 

64 . 0 

69 . 

7  i  .  0 

74.0 

7.  CRAQS 

13,0 

18.0 

L  8 . 0 

1 — 

CO 

3 

o 

18.0 

18.0 

25 . 0 

25.0 

25 . 0 

25.0 

25,0 

W  IN  CRACK 

18.0 

18.0 

18.0 

18.0 

18.0 

18.0 

25 . 0 

25. 0 

25.0 

25 , 0 

25 , 0 

BOUNDARY  1 

18.0 

18.0 

18.0 

18.0 

18.0 

18.0 

25 . 0 

25.0 

23 .  <:« 

25 . 0 

25.  0 

BOUNDARY  2 

25 . 0 

25 . 0 

25.0 

25.0 

25.0 

25 . 0 

25.0 

'•"kcr, 

25.  0 

25 . 0 

25.0 

BOUNDARY  3 

35 . 0 

35. 0 

25.  0 

25.0 

25 . 0 

23 . 0 

25,0 

25 . 0 

25 . 0 

22 0 

35,0 

BOUNDARY  4 

35.  0 

35 . 0 

20 . 0 

30.  0 

35.0 

35 . 0 

35.0 

35  0 

70.  0 

28.0 

35 0 

BOUNDARY  5 

75 . 0 

35.  0 

33.0 

■•‘U  »  0 

35.0 

35.0 

35 , 0 

35 . 0 

“’■tr 

OU  ..  '  * 

"■5  „  n 

32 ..  0 

REINFORCING  FORCE  DATA  AT  2  LEVEL (s) 


Y“  22.00 

NO.  OF 

X 

FORCE 

77 . 7 

.  0 

34 . 0 

5000 . 0 

“7  A 

/  . 

5000. 0 

.  0 

.  0 

Y--  23.00 

NO  OF 

Y 

FORCE 

'L2 . 0 

.  0 

39 . 0 

5000  •  '0 

7 , 0 

5000 . 0 

,  0 

,  o 

.  PROPERTIES 

COHESION 


|MI_  „  i  i 

i  00 . 0 


FRICTION  AND! 
70 . 0 
40.,  ■ 1 


DCHSIT 
i  15. ' 

l  r  •  >  - 

1  05 , 0 


1 000 .  0  40 . 0  1 20 . 0 
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NUMBER 

TANGENT 

RADIUS 

(X)  CENTER 

( Y>  CENTER 

F5 (BISHOP;  RS 

■.QMS)  c!FS  (38 

1 

35.  ^ 

13. 9 

37 . 0 

22 . 0 

CENTER 

BELOW 

INTERPOLATED  CP 

37.8 

l  5 . 8 

33 . 0 

22.0 

CENTER 

PELOW 

INTERPOLATED  CF 

35 . 7 

17.7 

37 . 0 

1 0 . 0 

.  863 

,  3**1  'j  „  40  * 

4 

35 . 

13.0 

4 1 . 0 

27  „  0 

CENTER 

BELOW 

INTERPOLATED  CP 

5 

36.  3 

10.3 

37 .  O 

26 . 0 

CENTER 

BELOW 

RIGHT  INTERSECT 

6 

36. 4 

18.  4 

35 . 0 

J  8 . 0 

..  “'81 

.  764  ..  377 

7 

35.  6 

19.6 

37.0 

1 6 ..  0 

.  032 

.911  .471 

8 

■-■5 . 3 

17.3 

39 . 0 

18.0 

.  759 

. 929  . 440 

9 

35.8 

15.3 

37 . 0 

20 .  'I1 

CENTER 

BELOW 

INTERPOLATED  CP, 

10 

37.2 

19.2 

33 . 0 

18.0 

.726 

.711  . 35* 

11 

36.  2 

20 . 2 

35 . 0 

1 6 . 0 

.843 

.333  .137 

i.2 

36.  6 

16.6 

35 . 0 

2C* .  O 

CENTER 
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INTERPOLATED  OR; 

13 

38.2 

20.  2 

31.0 

18.0 

2.729 

•T 

2  570  ,332 

14 

37 . 0 
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33 . 0 

16.0 

.  793 

.  730  ..  4 1  4 

15 

37 . 5 

17.5 

33 . 0 

20 . 0 
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BELOW 
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t  6 
17 
J8 
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36.  6 
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BISHOP  MODIFIED  AND /OP  CRD I MARY  METHOD  OF  SLIT 


permafrost  depth  7  feet,  10  feet  from  centerline 


CONTROL  DATA 
NUMBER  OF 


SPECIFIED  CENTERS 
NUMBER  OF  DEPTH  LIMITING  TANGENTS 
NUMBER  OF  VERTICAL  SECTIONS 
NUMBER  OF  SOIL  LAYER  BOUNDARIES 
NUMBER  OF  PORE  PRESSURE  LINES 

POINTS  DEFINING  COHESION  PROFILE 


It 


J 

0 

o 


NUMBER  OF 

SEISMIC  COEFFICIENT  S1,S2  =  .00,  .00 

UNIT  WEIGHT  OF  WATER  -  62.40 

SEARCH  IS  BASED  ON  BISHOP  MODIFIED  METHOD 
SEARCH  STARTS  AT  CENTER  (  37.0,  22.0)  WITH  FINAL  GRID  OF  2.0 

ALL  CIRCLES  PASS  THROUGH  THE  POINT  (  42.0,  35.0) 


GEOMETRY 


REINFORCING  FORCE  DATA  fiT  3  LEVEL (e) 


Y« 


20 . 00 
Y 

34 .  S 
3 1 . 0 
3 . 0 
,  0 


NO.  OF  FORCE  POINTS 
FORCE 

.  0 

1 666 .  (I1 
1666.0 


=  4 


Vir: 


22  .  00 


NO.  OF  FORCE  POINTS 
FORCE 


-=  A 


-■/  .  / 

r*  ^  *.) 


50^*0  -  0 


5000. O 


25 .  OO 

v 

42.  0 
39 ..  0 

2  .  o 


NO.  OF  FORCF  POINTS 
FOFCE 


5000 . 0 
5000.  0 

.  0 


a 


a 
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SECTIONS 

.  0 

1 0 .  0 

1 3 . 0 

15.0 

20 . 0 

72.0 

42.0 

64 . 0 

69.0 

71.0 

74. 0 

T.  GRACH’S 

18.0 

18.0 

18.0 

18.0 

18.0 

18.0 

25.0 

25 0 

23. 0 

23.0 

25 0 

W  IN  CRACK 

18.0 

18.0 

18.  0 

18.0 

18.0 

18.0 

25.0 

25 . 0 

25.,  0 

25.  0 

25..  0 

BOUNDARY 

1 

18.0 

18.0 

18.0 

18.0 

1 8 . 0 

18.0 

25.0 

25 . 0 

25.  0 

25 . 0 

25.0 

BOUNDARY 

2 

25 . 0 

25.0 

25.0 

25.  0 

25 . 0 

25.0 

25.  0 

25.0 

25 . 0 

25 . 0 

25.0 

BOUNDARY 

— r 

O 

35.0 

35.  0 

25.0 

25.0 

25.0 

25.0 

25.0 

25 . 0 

UmSmi  X 

n”  o 

-T-W  V  V 

35 0 

BOUNDARY 

4 

35.  0 

35.0 

28 . 0 

30. 0 

35 . 0 

35 . 0 

35 . 0 

35 . 0 

70  -  0 

29  4  0 

35 . 0 

BOUNDARY 

5 

35 . 0 

35.  0 

35.0 

35 . 0 

35.0 

35 . 0 

35 . 0 

35 . 0 

33 0 

35 . 0 

35.0 

► 
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SOIL  PROPERTIES 


LAYER 

1 


4 


COHESION 

.  0 

l  000 .  f'1 
100.0 
1000  .  0 


HR I CT I ON  ANGLE 
35.0 
40 .  0 


-1 0  .  0 


DENSITY 
135.0 
120.  :* 

1 05 . 0 
120.  0 


BISHOP  MODIFIED  AND/OR  ORDINARY  METHOD  OF  SLICES 


prr.na-f ro'it  depth  7  feet,  10  -feet  from  centerline 
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R  TANGENT 

RADIUS 

(X)  CENTER 

< Y  >  CENTER 
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1 

35.9 
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22.  0 
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n 
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33 . 0 

22 . 0 

CENTER 
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*T 
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37 . 0 
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.  887 
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37 . 0 

26 . 0 
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6 
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18.0 
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7 
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16.6 
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APPENDIX  C 


PROCUREMENT  OF  SSTIPN  PROGRAM 

The  procurement  of  computer  software  at  Georgia  Tech  must  follow 
Electronic  Data  Processing  Equipment  Request  Procedures  (EDP)  under  an 
executive  order  from  the  Governor.  The  procurement  of  SSTIPN  is 
summarized  as  follows: 

1.  EPD  request  number  11 -EC-86  was  submitted  for  approval 
(enclosed) . 

2.  Approval  letter  for  EDP  request  number  11 -EC-86  (enclosed). 

3.  Georgia  Institute  of  Technology  requisition  form  for  SSTIPN 
(enclosed) . 

4.  Letter  of  request  for  purchase  of  SSTIPN  to  Professor  Duncan 
(enclosed) . 


5.  Acceptance  letter  Professor  Duncan  (enclosed. 


Best  Available  Copy 


OFFICE  OF  ASSISTANT  VICE  CHANCELLOR 
FOR  COMPUTING  SYSTEMS 


BOARD  OF  REGENTS  OF  THE  UNIVERSITY  SYSTEM  OF  GEORGIA 

PRESIDENT'S  OFFICE  RUILOING 
ATHENS.  GEORGIA  30*02 


January  31,  1986 


Dr.  Joseph  M.  Pettit 
President 

Georgia  Institute  of  Technology 
225  North  Avenue,  N.W. 

Atlanta,  GA  30332 

Dear  President  Pettit: 

Approval  is  given  to  Georgia  Institute  of  Technology  to  acquire  the 
data  processing  equipment  referenced  in  request  number(s): 


Yours  very  truly, 

•  -  r> 


{sCsy'Y’Jt’iL, 


.  -v  to-yV's. 


James  L.  CArmon 
Assistant  Vice  Chancellor 
for  Computing  Systems 


JLC/jw 


cc:  Mr.  John  Gehl 


CEORGIA  INSTITUTE  OF  TECHNOLOGY 


SCHOOL  OF  CIVIL  ENGINEERING 


-2-2-  /£4> 
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Signacure  of  person  requcscing_ 


■y-'./  lLcX-^ 
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School  of  Civil  Engineering 
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Best  Available  Copy 


COLLEGE  OF  ENGINEERING 


VIRGINIA  POLYTECHNIC  INSTITUTE  AND  STATE  UNIVERSITY 

* 

Blacksburg,  Virginia  24061 

DEPARTMENT  OF  CIVIL  ENGINEERING  (703)  961-6635 

February  13,  1986 


INVOICE 

Mr.  George  Papaioanou 
c/o  Dr.  Neil  D.  Williams 
Georgia  Institute  of  Technology 
Department  of  Civil  Engineering 
Atlanta,  GA  30332 

REF:  Requisition  No.  200-6-25147 

For  computer  program  SSTIPN  (Soil-Structure  Interaction  Program) 
sent  via  BITNET  on  February  13,  1986. 

SSTIPN  $200.00 

Payable  upon  receipt. 

Mail  payment  to: 

Professor  J.  M.  Duncan 
Department  of  Civil  Engineering 
104  Patton  Hall 
Virginia  Tech 
Blacksburg,  VA  24061 


Please  make  check  payable  to:  The  Virginia  Tech  Foundation,  Inc. 


APPENDIX  D 


BASIC  PLOTTING  PROGRAM 

The  BASIC  program  for  plotting  FEADAM84  and  SSTIPN  reduced 
output  was  written  by  Robert  L.  Roglin  (21).  The  program  was  written 
for  use  on  an  IBM  PC  with  enhanced  computer  graphics.  Output  data  from 
FEADAM84  and  SSTIPN  was  reduced  for  use  in  the  plotting  program.  Four 
data  files  were  reduced  from  the  output:  Nodal  point  input  data,  four 
node  solid  element  data.  Nodal  point  displacement  data,  and  four  node 
element  stress  data.  Three  dimensional  plots  of  stresses  were  not 
useful  in  evaluating  stress  distributions  within  the  embankment  and  are 
not  included.  Program  source  listing  is  enclosed. 


LIME 

"'EM 

REM 


(o.o)  --  •:  4-ot  „  ■)) 


[  ME’  <  0  ,  >.') 


'  :>,  if.- 


X4  = 
LINE 
NEXT 
REM 
REM 
REM  I 
REM 
REM  - 


1 t  L n de f 01  mad  Goo-ner.r> 


=  i  TO  MEL 

HOC  ALE*  '<  ( I  L  (  EL  ( I  >  >  ;  ;  Vi 
HOC  ALE  ( '<  ill:'.  EL  ilu;  r,  V  2 
HSCALE* X  *  E 3 (EL  1 1 > ) ) j  Y3 
HSCALE* X  t 1 4 ( EL '!>>>:  Y4 
-  (XI, Yi) 


Plot:  De-f  armed  Geometry 


VSCAi-E  *  Y  < 

/ji.hiulli  l  i 

VSCALE  Jt  Y  t 
V3CALEYY • 


i  l  1  C.L  L  1  I  / 

12;  IL  - 1  >  • ; 
1 3 { EL ( I >  )  > 
[ 4 ( EL  v I ) ) > 


_ L  i  .  i _ 

LIME 

LINE 


37  V  r\L.:  i 

900  FOR 
?  1 0  X 1 


FOR  I  =  1  TO  NJT 

X 1  =  mSCALE t X (1 1 (EL  il ) ) )  *  DHSCALESDX (LC . 
X2  =  HSCAuE  KX  (  J.2  ■  EL  c I )})  +  DHSCALEfrD*  (LC, 
X3  =  HSCALE* X  (  E3  (EL  ( I  >  >  )  +  OHSCALESDX  ■'  LC. 
X4  HGCALE*X  ' 14 (EL  (I >  >  >  *•  DHSCALESDX  (LC, 

Yt  =  VSCALE*  Y  ( I  l  (EL  ( I  >  )  )  +  DVSCALEiDY  .’LC, 
Y2  =  VSCALESY(I2'EL(I) >>  +  DVSCALE*DV (LC, 
Y3  =  VSCALE *Y (13 (EL ( I > ) )  +  DVSCALEIDY (LC, 
Y4  =  VSCALE* Y  f 14 (EL ( I )  )  )  +  DVSCALEYDY (LC, 


(EL (I  1 > / 
:  ■  EL  ( I  :• ) ) 
(EL  (I)  ;  ■ 
• . EL (I  •  )  > 


990 

LIME  (XI 

■ ,  Y 1 ) 

- 

■.  C  00 

LINE  - 

( X  3 , 

Y3 ) 

10  10 

LINE  - 

1X4, 

Y4) 

i  020 

LINE  - 

( X  i , 

Yl) 

i  OZ0 

NEXT  I 

l04i() 

REM 

1 050 

REM 

C'  RhM  Stress  Con t our  F 1  ot t i r. g 
0  F'Ei’i 

3  REM  - 

0  REM 
0  REM 

0  REM  Eevelop  Coordinate  Transformation  Ei 
2  REM 

I)  FEM  - 

0  REM 

0  INPUT  "SPIM  ANGLE,  TIP'  ANGLE:  " :  IN,  T 1 
0  CCPIM  =  CCS«3PIN*3.  H 155/1  F.C) -•  ES^'IN  = 

0  CTI*  =  CCS  ( T I c'  *  3 .  1 4 1 5  9 1 30  '■  :  ST  (P  ■=•  3,.' 

0  PRINT  "cap  in  =  "  .CSPIN,"  ssn-r*  -  “.TIP 
0  -PINT  "ctio  =  ",cr:p.  -  3t •  _•  -  ■ .  3 " I P 

INPUT  "STRESS.  LOAC:  " :  CTREE"  . .  Z 


2  : 

rTr  r~ - 
.*>  I  »’  ll  L/O 


•'•rr, 

>  >  i _ — 

1  THEN 

2  T  :jEM 
~  ~FEN 


MEN  2  L-5  --  3* 
EL  I  >  MC 
EL.  ’  .  )  •  r.-’i  ;  r- 


T  v  •  l'i  :-L 


‘  ME  X' 


REM 

"l  -“-il 


t  L tress  Cont  in.  s 


|  1330 

[  1340 

L  1 350 


it*  i,-*  vw  v*  v>i  v^,vvr»^.uv^-u«.w?vi)VT^ri.V7Aa^CT«v%*JWii-VY’c 


CTXTnrrrj  »v.*.v 


1 3V0 
!  380 
1 390 
1 400 
1410 
1 4  20 
1 4  3<  > 
1440 
1450 
14o0 
l  430 
.1  400 
14  90 
1500 
1510 
1520 
1 530 
1540 
1 350 
1560 
L  570 
1 580 
1590 
J  6  00 
1 6 .1 0 
1620 
1 630 
1640 
1650 


Plot  Globa,.  A; 


=  LOOiCSPIN:  Y3X  =  -1  OO'^SPIN-SG'^IP;  -.[WE  ( 
=  0:  7SY  -  100*01  IP:  LINE  <0,V)  --  ( XS7 ,  YSY 
*  iOO#3SPE!M:  YGZ  ==•  1  .'04CSP  TNKSTI F:  LINE  <0 
:  -•  1  TO  MEL 

hSCALP*X  ■;  1 1  (EL  (  J  >  >  }  5  V  i.  -=  VSCALESY  1 1 1  (EL  ( 
YltCTIP  -i-  Xl*3TIPlc(-l  4S0PIN:  XI  =  XL*CSF 
HSCALE  *  X  ( 1 2  ( EL  ( I  >  )  )  s  Y2  ==  VSCALE* Y  ( 1 2  ( EL  < 
i'ZUCTIP  +  X2*STIPf  (~t )  CSSPIM:  X2  =  X2*CSP 
HSCALE «  < I  3 (EL ( I )  )  >  :  Y3  =  VSCALE * V (1 3 ( EL  v 
Y3KCTJP  +  X3*STIP#(-l)*SSPINs  X3  =  X3SCSP 
HSCALE  It  X  (1 4  ( EL  (I )  )  )  ;  Y  4  =  VSCALE  K  Y  (1 4  ( EL  ( 
Y 4 * CT  IP  +  X4K ST  IPX  ( ~  1 )  % SSP IN:  ,14  =  X 4 4 CSP 


XSX  „  Y* 


\  v2  i  i  •  J  *„  , 


IN 
I)  >  ) 
IN:  L 

n )  > 

I N ;  L 
C)  )  * 
IN:  L 


<xi,  1 1 )  -  \ ; 


LINE  -  (  X  1  j,  Y 1 ) 
NEXT  I 

FOR  K  »  L  TO  NLR 
ia  ^  ( NLR  -i-  1 1  - 


-  ELS  (Li)  TO  (ELF  (1  1 )  -  1 
ELF  (I'M  )  -  J:  IE  =  ELF  (H  I > 


J  -  1 


LINE  (XS(IB)  VS  (IB)  )  -  <  XS  (IE)  ,  YS  ( IE)  J 
IF  HI  -•■=  1  THEN  GOTO  1640 
IE  «  ELF (K1  -  1)  -  J 

LINE  (XSlIB)  ,  YSKB)  >  -  (XS  ( IE)  ,  YS  (  TEM 
J  =  J  +  1 
NEXT  T 


1660  IF  K 1  -  1  THEN  GOTO  1700 
1670  IB  =  ELF (HI)  -  J:  IE  =  ELF (LI  -  l>  -  J 
L<sSO  LIME  (XS(IB)  ,  YSdB."  -  (  *S  i  IE)  ,  VS  ( IE)  / 
L690  NEXT  H 
L700  REN 
L  1710  STOP 


APPENDIX  E 


SSTIPN  OUTPUT  DATA 

Detailed  results  of  the  SSTIPN  analyses  are  contained  herein. 
Output  data  includes  nodal  point  displacements,  soil  element  stress  and 
strain  data,  and  internal  forces  of  structural  elements.  Analyses 
include  output  for  all  five  different  geotextile  strengths.  Each  of 
the  five  geotextile  strengths  include  three  different  conditions  as 


stated  in  Section  3.2. 
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UUMIfl" BF” BIFF .  HAM  HZmiALl-'- 

NUMOKR  OF  NOOAL  LINKS  . 

NUMOKR  OF  INTKRFACK  KLRMKNTS . 

NO  OF  INTKRFACK  RLK  IN  FRKRXIST  FART 

NUMSRR  OF  INTKRFACK  MATKRIALS . 

TOTAL  NUMOKR  OF  SOIL  KLRMKNTS . 7 

NUMOKR  OF  OIFF  SOIL  MATKRIALS 


NUMBIft  OF  IIIMIMTI  I»  “FUUMOATIBH  ” 

NUMKKR  OF  NOOKS  IN  FOUNDATION . 

NUMOKR  OF  FRRRX 1ST  INC  KLRMKNTS -  *  •  “ 
NUMKKR  OF  FRRKXISTINC  NOOKS . 

NUMKKR  OF  LOAD  CASKS . 
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101 

41 

49 

1  1  1 

314  C 

4  0 

411 

-  211 

•If 
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SO 
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12 

13 
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311.1 
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211.1 
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4.0 

4.0 
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.411 

.411 
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•  .221 

-  073 
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.310 

.211 


2.110 

2.711 

on 
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-  .211 
-.221 


.211 

-  .  Ill 

013 

v/aar 

3.107 

4.433 


2.111 

2.711 

074 
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1.411 

2.227 


•  .333 
►  .221 

•  077 

-YVfU” 

-1.401 

-2.221 


3.201 
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3.112 

4.411 
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12 
13 
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It 

11 


I 


17 

1 
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O 
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013 
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17 
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— rm — 

.077 

- TXT. 

3. >71 

r.  (T  r 

3.431 

'  -.IH 
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It 

11 

1213.2 

711.1 

•  17.4 

.  133 

-  .312 

2.713 

•  .  341 

2  713 

-  .31* 
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tt 

•  0 

411.3 

411  .  1 

111  .  1 

.211 

-  .  102 

.014 

.027 

.011 
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•  1 

11 

.  1 

4.0 

O 

.411 

.171 

-.121 

2.474 

4.310 

-1.330 

2.710 

•  1 

•  2 

.  1 

4.0 

.0 

.411 

.344 

-.231 

1  .207 

.723 

-  .117 

1.340 

12 

•  3 

211 .7 

430.0 

14. 1 

.311 

.221 

-  .  171 

.311 

.213 

-  .221 

.sot 

13 

•  4 

.  1 

4.0 

.O 

.  411 

.321 

1.121 

11.111 

14.211 

•4.017 

11.341 

•  4 

•  1 

.  1 

4.0 

.  o 
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-1. 444 

23.41« 

•23.141 

27.111 

•  10.122 

31.201 

•  1 

•  1 

.  1 

4.0 

.  0 

.411 

-2.410 

11.211 

♦14.312 

14.240 

-41.101 

•  1.040 

11 
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1 

4  O 

o 

411 

•  114 

•3  HO 

•17  211 

30  SOI 

-27  111 

11  214 

•  7 

“  tt 

.  1 

4  O 

.  6 

.111 

.  (31 

‘ - .6(5 

.  STM 

'  'IT!" 

• .  6  'i  ( 

.311 

Cl 

•  1 

212.9 

423.1 

71.4 

.411 

.234 

-.211 

.213 

.211 

•  .243 

.SCI 

•  1 

70 

.  1 

4.0 

.O 

.411 

.111 

10.137 

•1 . *71 

10.172 

.123 

1.141 

70 

71 

1 

4  O 

O 

411 

-12  207 

21.137 

-13  011 

27  10O 

-13.271 

41  171 

71 

. 

. v . 

. r/6 . 

-V4»«' 

•••-37343' 

-V77S3V 

'"ii  :itv 

347464 

“  »v.  nr 

16.  (i’i’6 

73 

.  1 

4 .0 

.0 

.  411 

11.731 

-10.211 

*21  117 

21 .ill 

•11.411 

37.411 

73 

74 

.  1 

4.0 

.0 

.411 

.210 

-  .  113 

-7.411 

3.711 

•3. Ill 

7.411 

74 

71 

331  1 

431  1 

123  S 

311 

OOI 

OlO 

Oil 

01  1 

001 

012 

71 

— Tt 

"■  TT 

4  .  O 

- "5 - 

.411 

•  n  3ii 

11.443 

4  n  i 

II. IH 

-11.421 

30  020 

71 

77 

.  1 

4.0 

.  0 

.  411 

• tO. 131 

11.311 

-1.344 

11  747 

-  10.413 

21  230 

77 

71 

.  1 

4 . 0 

.0 

.411 

20.073 

-14.171 

-  .0.112 

20.171 

-It. C77 

31.313 

71 

i 


FOUR  NOOE  SOL  10  ELEMENTS  •  STRESSES 


"TCI - 

HR - 

— rm — 

TAU-HV - 

tie-i 

— rrr-s 

tAu-ma4 

- THTTA 

nciyrro — 

— trar 

ELE 

1 

.313 

.114 

.ooi 

.114 

.213 

.  121 

I  .  121 

1  .111 

on 

1 

2 

342 

417 

040 

103 

23f 

133 

1  141 

2  131 

on 

2 

. 3 . 

. Til* . 

. '.•3T3 . 

. ."(34 . 

. :t*4 . 

. .'Y3C  ’  ’ 

'"*627311" 

.  a. 36*' . 

.’613* 

. 3 

4 

.313 

■  14 

.  001 

.  114 

.313 

OH 

1  .  92 1 

1  447 

1  717 

4 

1 

.371 

ISO 

.on 

.111 

.370 

.Ot  1 

4  .  Ill 

1  411 

1.112 

1 

• 

413 

•  It 

013 

•  It 

411 

on 

3  112 

1  424 

1  III 

• 

7 

- TT5 - 

.7  12 

.  007 

""  "7TV5 

.  ITi" ' 

TT53 

i  '134 

( .313 

3.oo3 

7 

1 

.333 

•  S23 

-  .  002 

.  123 

.  333 

.011 

•  121 

1 .172 

1  101 

• 

• 

.310 

.117 

.023 

.111 

.311 

.311 

2  037 

2.771 

in 

• 

10 

401 

1  447 

021 

1  441 

401 

122 

1  112 

..  .V4  ;m 

3  171 

210 
*  -311 

10 

. n . 

.  446 

'(".■436 . 

. fit  "" 

"  i  ."414  ‘ 

. 4  VS" 

. V«Y* ' 

. 3 . *** 

12 

311 

.  110 

-  100 

.121 

.  312 

.  213 

•  10.311 

2.110 

114 

12 

13 

.  Ill 

.301 

.ooi 

.  301 

.  Ill 

on 

3 . 101 

1 .130 

042 

13 

14 

134 

211 

017 

213 

on 

071 

21  413 

2  121 

013 

14 

(I  ' 

- 7TT3 - 

.  174 
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TTTT - 

:  too 

.031 

io  TVo 

1  .  31  (  ' 

714 

II 

11 

.271 

.211 

.031 

321 

.  243 

.031 

40.401 

1 .311 
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It 

17 

.334 

.417 

017 

420 

.  331 

.  041 

ii  no 

1  .  270 

114 

17 

11 

312 

411 

oto 

411 

311 

014 

1  107 

1  301 

. f  331" 

1  017 

11 

•"•(*  "  - 

■".‘464"” 

6(3 

. ."411  ' 

. iti"" 

.  663  ’ 

'  3  166 

20 

412 

.101 

.  009 

•os 

.  41  1 

.014 

4  727 

1  231 

1  113 

30 

21 

.310 

417 

.  OOI 

4  17 

.310 

.014 

2  711 

1  341 

1  071 

21 

22 

331 

1  111 

OlO 

1  170 

333 

•  11 

2  S7f 

1  001 
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22 

'  33  ' 

.  323 

— rrm — 

- -  I’M - 

— rrm — 

.  211 

- .471“” 

-io.iW 
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301 

33 

24 

.319 

.  ill 
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.  117 

.  217 
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•If.  lot 

2.421 
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24 

21 
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on 

•0.711 

S.401 

.071 

21 

21 
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-  007 
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•1  017 

"4366’.  (ir  * 
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”•-"(  666' 
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-■nr  • 
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.  666 
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. .611  " 

.  006 
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•••-»*•  iti 

21 

.117 

.  034 
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.117 

.  034 

.042 

•n  .no 
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•  34 

21 

21 

.  112 
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.  214 
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1 .170 
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21 

30 
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•  002 
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m 
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*  441 
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30 

31 
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31 

32 
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.  210 

-  013 

.  21 1 

no 

.  Of  t 
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32 

33 

.414 

.  101 
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.  414 

.017 
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1.231 

1.131 

33 

34 

.  ••• 
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111 

in 
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on 
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34 
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31 

31 

.017 

.304 

• .  001 
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.  017 
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4.  S44 

on 

31 

37 

.377 
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•in 

.  411 

.  Oil 

220 
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37 

31 
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on 
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on 
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•  770 

2  111 

31 

31 
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.  ok  i 

1  lo 

- 33.113 

1  '4  IK 

■3.771" 

31 

40 

•  .202 

.227 

•  121 

211 

•  234 

241 

-14  713 

*  1  107 

•  1  OOO 

40 

4  1 

*  217 

.121 

•  .331 

.  114 

•  .  312 

111 

-20.313 

•  1  710 

•  1  OOO 

4  1 

43 

•  .714 

.333 

.111 

341 

•  .  7*7 

171 

1.117 

•  .433 

-  1  .OOO 
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r  1 

I 

w  i 
■  4 

•  i 
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& 
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SIT  I  PM :  t  LAYCR  OP  CCOTCXT l L E ,  P«2000  PPt  ,  T-0 

TOTAL  NUMBCR  OP  NOOCS .  SS 

NUMBCR  OP  BAR  CLKMKNTS .  t 

NUMBCR  OP  OIPP.  BAR  MATIRIALS . .  1 

NUMBCR  OP  BCAM  CLCMCNTS .  O 


jiri  i  I  BIJI'HJJW  1  |.|f  ^ 


NUMBCR  OP  NODAL  LINKS  .  < 

NUMBCR  OP  I NTCRPACC  CL CMC NTS . 

NO  OP  INTKRPACC  CLI  -  IN  PRBBXIST  PART 

NUMBCR  OP  J NTCRPACC  MATBRIALS . 

TOTAL  NUMBCR  OP  SOIL  CLCMCNTS .  T 

NUMBCR  OP  PIPP  BOIL  MATCRIALS . 


1*111  1  IN  411  1  R.'  l'M-Mir 


NUMBCR  OP  NOOCS  IN  POUNOATtON . S2 

NUMBCR  OP  PRCCX ISTINC  CLCMCNTS .  O 

NUMBCR  OP  PRCCXISTINC  NOOCS .  O 

WU»rf||fk*'dP"eO«DYKU<iTttf«"lAVOB--vvv:-;-:-vv-"*4 

NUMBCR  OP  LOAD  CASCS .  1 

CAL  INC  PACTOR . . . 1  OOOOO _ 


ATMOSPHCR I C  PRCSSURC  •••  1 . OSSOO 

UNIT  WCI6MT  OP  WATCR -  0312C 


m  w  RiiiXBin.au 


i 
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PUT  ON  LAVCR  NO.  2 


NODAL  POINT  INPUT  OATA 


NODE  NODAL  POINT  COORDINATES 
NUMtEfc  X  *  0N0  Y  *  0 AO 


10.000 

4.000 

20.000 

4.000 

30  000 

4  OOO 

31.666  4 .666 

42.000 

4.000 

•0.000 

4.0CO 

•  4.000 

4  OOO 

44T6oo 

47663 

•  •.000 

4.000 

74.000 

4.000 

•2  000 

4  OOO 

- 17 35? 

.ooo 

7.000 

10.000 

7.000 

20.000 

7  OOO 

74.000 

t 1 .*00 

•2.000 

1  t  .too 

•0.000 

1  i  .too 

1  *2.000 
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t  1  soo 

i6  At. 6b 6 
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»  SO. ooo 

1  1  .*oo 
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42.000 

1 3 . OOO 

SO . OOO 

13.000 

•4.000 
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••.000 
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••.OOO 

13.000 

74  OOO 

13  OOO 
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14.000 
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ts.ooo 

4  OOO 
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4.666 

\ 4 . 666 

• .  ooo 

IS .ooo 

4.000 

IS .000 
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•  4.000 
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IS. ooo 
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•  •  ooo 

1  • .  ooo 
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14.666 
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10 . ooo 

IS. ooo 
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17  OOO 
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0 

o 

o 

4 

•2.000 

17.000 

0 

0 

o 

t 1 

SO  OOO 

17  OOO 

1 

o 

1 

»* 
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STRUCTURAL  ELEMENTS  *  LINEAR  ELASTIC 


YOUNG'S  MOOUl 

••nm  •■•uwrrwr  "'wimiir 


•OOO . JO 
»0  06 


•ULK  MODULUS 

Arid  **  ttumirr”  WR6<r*MV 


1100.00 
20 . 00 


STRENGTH  PARAMETERS 

**if . m  ”  ‘6R«r* 


KdO'A"MbbMI'“i‘6wYb"tL'fcrtKWY 'tSAt'A . 

SLIT  CONNECTED  KOOES  MATL  ELEMENT  CENTER  COORDINATES 
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2 

IS 
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2 

2 
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20 
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20 
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n 

13 
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IS 

24 

27 

14 

IS 

11 

2t 

20 

IS 

It 

17 

30 

2* 
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2.000 
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2.000 
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iY. bdd 
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41.000 

2.000 

S2.00C 
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S  .  soo 
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fl— ~  T1 1  T 1 
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21 

22 

3f 
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■ 
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22 

23 

3S 
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23 

24 
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23 

24 

2S 
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2S 
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27 

24 
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r 
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3S 
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34 

34 

34 
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S2 

so 

SI 

3 

3 
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37 

43 
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S3 

S3 

? 
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ji 
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SS 

S4 

1 
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40 
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47 

47 

4  A 
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SS 
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1 
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St 

so 
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70 
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44 . OOO 
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y 
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53 
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74 

*7 
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75 
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so 

SI 
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41 . bob 
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1*2 -  i 
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44 
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i 
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7S 

SS 
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SS 

7t 

SO 

St 
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• 
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■1 

71 

si 

42 

SS 
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72 

42 

S3 

ss 

SS  1 
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73 

S3 
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7S 

AS 

S  7 

42 
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S7 

SS 

43 
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1 S . SOO 
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77 

SS 

44 
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4o 

”  44  * 
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*% 

f. 
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CONSTRUCTION  ER  INFORMATION 
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•\«*  R  A 


ISJCTJTRy  r, kl*.k:£Jf35S214"*  *7n? V*  CT%7.  WTTW^^^rwirW.^'* . G"VrtTM  2.7 - 

LOAD  CASE  t  1  ITERATION  *  2 


NF 

DELTA- X 

OELTA-Y 

OELTA-22 

X-OISP 

Y-DISF 

22  - ROTAT 

TOTAL 

NP 

\ 

.  0066 

7 0066” 

.  0606 

,  0000 

aaao 

.  6606 

66  J« 

\ 

2 

.0000 

.  0000 

.0000 

.  0000 

.  0000 

.0000 

.0000 

2 

2 

.0000 

.  0000 

.  0000 

0000 

.0000 

.0000 

.6100 

3 

4 

0000 

0000 

0000 

0000 

0000 

0000 

.  0000 

4 

...  | .. 

. '6666'" 

'.'6600 

“*  .0666 '" 

;'o666'" 

"  ".‘6'6'66  ■ 

. .'6660 

. '.  6666 

.  .  g. 

• 

.0000 

.  0000 

.  0000 

.  0000 

.  0000 

.0000 

.0000 

8 

7 

.0000 

.  0000 

.  0000 

.  0000 

.  0000 

.0000 

.0000 

7 

t 

0000 

0000 

0000 

0000 

0000 

0000 

RHmai 

• 

i - 

.6606 

- 5756 - 

.  6666 

■t-i-m 

.0666 

■f-M-I-K 

■I'H-r-a 

r 

10 

.0000 

.  0000 

.  0000 

0000 

.  0000 

.  0000 

.  0000 

10 

1 1 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

.  0000 

.0000 

1 1 

12 

.  0020 

0000 

0000 

0000 

0000 

0000 

0090 

12 

■'i  J  ' 

.  6006 

6660 

.  0060 

.  6666 

,  6666 

.  6060 

.  0066 

‘  '1'3 

14 

.0000 

0000 

.<*000 

.  0000 

.  0000 

.  0000 

.0000 

14 

It 

.0000 

0000 

.0000 

-  .  0002 

.  0000 
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MATL 

UNIT  WT 

CONSTANT 

EXPONENT 
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